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REMARKS 

Dependent claim 50 has been amended to depend from claims 42-44 and 47-49, since 
claims 45 and 46 were cancelled. Claim 59 has been cancelled without prejudice or disclaimer. 
No new matter has been added. It is believed that the entry of the amendments to claim 50 and 
the cancellation of claim 59 would not require a new search and examination to be performed, 
and thus, it is respectfully requested that this amendment be entered. 

Claims 33-35, 38-44, 47-50, 56-58, 61-63, 70 and 72-77 are currently pending for 
examination. 

Rejections under 35 U.S.C. $101/$! IZ HI 
Claims 33-35, 38-44, 47-50, 56-59, 61-63, 70 and 72-77 have been rejected under 35 
U.S.C. §101, as the Patent Office states that the disclosed invention is inoperative for the full 
scope of the claims and therefore lacks utility. These claims have also been rejected under 35 
U.S.C. §112, T[l , as failing to comply with the enablement requirement, for similar reasons. The 
Patent Office asserts that Hirvonen, et al "Effect of Diffusion Potential, Osmosis and Ion- 
Exchange on Transdermal Drug Deliver, Theory and Experiments," J. Controlled Re I. 56:33-39 
(1998) ("Hirvonen") provides evidence of the operability and unenabling of the present invention 
as claimed. 

Initially, it should be noted that the claims as pending each recite an L-arginine or an L- 
arginine derivative that is applied to skin. The Applicant is not claiming that any possible 
substance can be delivered to the skin under the conditions specified in the claims, but rather, 
only the L-arginine or L-arginine derivatives as claimed. As evidence, the Applicant has 
provided evidence, both in the form of a declaration as well as a series of examples originally 
filed with the patent application, to show that L-arginine or L-arginine derivatives can be applied 
to the skin under such conditions. 

The Patent Office has not provided any evidence that illustrates that L-arginine or L- 
arginine derivative cannot be successfully applied to the skin that contradicts the Applicant's 
evidence. Instead, the Patent Office has pointed to Hirvonen as a reference that the Patent Office 
asserts that one or ordinary skill in the art, in reviewing the specification, would conclude that L- 
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arginine and L-arginine derivatives cannot be applied to the skin of a subject as directed by the 
specification, despite the fact that there are working examples in the specification that teach 
otherwise. Hirvonen nowhere discloses or suggests the use of L-arginine or L-arginine 
derivative, and thus its usefulness in this regards is limited, especially in comparison with the 
specific evidence in the application as filed to the contrary. 

In particular, Hirvonen teaches that sotalol cannot be delivered to the skin under certain 
conditions. The Patent Office then extrapolates, based on the sotalol data of Hirvonen, that one 
of ordinary skill in the art would then conclude that the invention as presently claimed is 
inoperable and non-enabled. However, Hirvonen is inconsistent with itself^ thus its validity as a 
suitable reference would be highly questionable to one of ordinary skill in the art. In particular, 
Hirvonen teaches in Table 1 that for Composition VI, theoretically, there should be relatively 
little transport, as Hirvonen predicts Edp = 1.15, where Edp is the enhancement factor of the 
increase in transport relative to controls where there is no osmotic difference. However, the 
value actually measured (E"dp) was 8.40 ± 2.74 (SE), which is not only statistically significantly 
greater than the theoretical value of 1 . 15 (by more than 2 SE units), this value was the largest 
enhancement factor actually measured by Hirvonen on Table 1 . Thus, not only the validity of 
Hirvonen' s theoretical treatment rendered suspect by its ovm data; Hirvonen' s data also actually 
shows that even for sotalol, changes in osmolarity can affect molecular transport, in contrast to 
the position of the Patent Office. 

Moreover, the formulas for sotalol and for salicylate, both discussed in Hirvonen, are not 
structurally similar to L-arginine or its derivatives. The formula for sotalol is shown below: 




H3C O 



HN 




and for salicylate: 
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while the structure of L-arginine is as follows: 



H 




The Patent Office has not asserted that any of these chemical formulas are structurally 
equivalent, nor has the Patent Office has not pointed to any references that teach this, or suggest 
that any results obtained for sotalol or salicylate can be applied to L-arginine or L-arginine 
derivatives. Given the uncertainties inherent in transdermal drug delivery, as evidenced by 
Hirvonen (among other references in this field of study), one of ordinary skill in the art, in 
examining the experimental data demonstrated in Hirvonen, would not believe this data and the 
experimental data within the instant application to be so totally inconsistent that this person 
would then conclude that all of the Applicant's data does not support that L-arginine can be 
delivered as claimed, and all of the instant claims as pending are therefore invalid for being 
inoperative for their full scope. 

Moreover, research into the scientific literature concerning Hirvonen reveals that 
Hirvonen itself is not a widely-accepted paper. In fact, to date. Applicant has been able to 
identify only three papers have cited Hirvonen, after a study using SciSearch. These three papers 
are Haldiya, et al, "Dermal Ulcers and Hypertension in Salt Workers," Current Sci,, 87:1 139 
(2004) ("Haldiya") (Appendix 1); Suhonen, et al^" Epidermal Cell Culture Model Derived from 
Rat Keratinocytes with Permeability Characteristics Comparable to Human Cadaver Skin," Eu. 
J. Pharm, ScL, 20:107 (2003) ("Suhonen") (Appendix 2); and Matuszak, et ai, "Thermodynamic 
Driving Force for Molecular Diffusion - Lattice Density Functional Theory Predictions," J, Non- 
Equilib. Thermodyn,, 31:355 (2006) ("Matuszak") (Appendix 3). 
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Suhonen is directed to a cell culture model, and only mentions Hirvonen in passing, and 
then only for the proposition that Hirvonen successfully delivered solatol across the skin (see 
Table 1, p. 109). Suhonen is concerned only with showing that solatol can likewise be delivered 
across cultured epidermal cells, rather than skin. 

Haldiya was a study of workers exposed to high salt concentrations (brine), who 
exhibited ulcers or hypertension. According to the paper, salt concentrations did positively 
affect the incidence of ulcers and hypertension. Such effects would not be observed if the salt 
concentrations had no effect on the skin. In this paper, Hirvonen was cited for the proposition 
that "When salt concentration in fluid surroxinding the skin is high, osmotic pressure and ion- 
exchange make a significant contribution to flux enhancement by diffusion potential" (p. 1 140, 
second column, first paragraph), contradicting the position of the Patent Office regarding 
Hirvonen. 

Matuszak appears to be the only paper that mentions the theoretical model of Hirvonen. 
However, Matuszak states that the diffusion potential gradient for diffusive flux in Hirvonen is 
only an approximation, and thus should not be treated as a universal linear law. See p. 359, 
paragraph 2. In other words, Matuszak states that Hirvonen' s theoretical equations have only 
limited applicability. In fact, the equations that were derived in the Matuszak reference appear to 
support the Applicant's findings that transport can occur through the skin under such conditions. 
Thus, a person of ordinary skill in the art that is aware of Hirvonen would also be aware of 
Matuszak, which corrects Hirvonen' s equations, and would accordingly not question, based on 
Hirvonen, whether agents for creating a hostile biophysical environment could be used to cause 
transport across the skin to occur. 

Accordingly, given that Hirvonen has been published now for more than 10 years, the 
relative lack of citations by others in the scientific community indicates that Hirvonen itself is 
not widely accepted by those of ordinary skill in the art. Papers that reference Hirvonen only 
refer to the experimental data or state that Hirvonen' s theoretical equations are not broadly 
applicable in all situations. Since Hirvonen's own data discredits those theoretical calculations, 
it should not be surprising that no other paper in the scientific literature has placed any 
significant weight on those calculations. Accordingly, as this evidence shows, one of ordinary 
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skill in the art would have no rational reason to believe any of the theoretical calculations 
discussed in Hirvonen in the face of contradictory experimental data (in both Hirvonen and in the 
instant application) and literature statements that do not support Hirvonen. Therefore, Hirvonen 
should not be relied upon to sustain this rejection under §101. 

Thus, for at least these reasons, it is respectfully requested that the rejection of these 
claims be withdrawn. 

Rejections under 35 U.S.C. $112. 112 

Claims 50 and 59 have been rejected under §112, 112. Claims 50 (in part) and 59 depend 
from cancelled claims. 

Claim 59 has been cancelled without prejudice or disclaimer. Claim 50 is a multiply- 
dependent claim, and the dependencies in claim 50 have been corrected. Thus, it is respectfully 
requested that the rejection of claims 50 and 59 be withdrawn. 

Declaration 

In rejecting the declaration as being insufficient to overcome the rejection of the claims, 
as discussed in the Office Action, the Patent Office appears to put greater weight on the 
theoretical treatment of Hirvonen than on the Applicant's data in both the application and in the 
declaration. For instance, the Patent Office states that Hirvonen teaches that different 
mechanism act to defeat enhancement of transdermal drug delivery at high and low electrolyte 
concentrations, and that Hirvonen provides a detailed theoretical discussion as to why in general 
enhancement of transdermal drug delivery by added electrolytes should fail and for this reason it 
must be assumed that in general such enhancement should fail. 

However, Hirvonen teaches that this is true only for sotalol and salicylate; Hirvonen 
offers no data or mechanisms for L-arginine and L-arginine derivatives, and Hirvonen's 
theoretical model fails even its own experimental data for predicting transport of sotalol and 
salicylate, as discussed above, to say nothing of L-arginine and L-arginine derivatives. As 
Hirvonen's own "detailed theoretical discussion" emphasized by the Patent Office is 
contradicted hy Hirvonen's own experimental data, Hirvonen's detailed theoretical discussion is 
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not useful for predicting the transport properties of L-arginine and L-arginine derivatives (or 
indeed, for even predicting the transport of sotalol or salicylate), and accordingly, Hirvonen's 
failed theoretical model cannot be used to contradict the actual experimental data obtained by 
the Applicant. 

Rejections under 35 U.S.C. gl02(bV$103(a) 
Claims 33 and 61 have been rejected under 35 U.S.C. §102(b) as anticipated by, or in the 
alternative, under 35 U.S.C. § 103(a), as being obvious over Ermen, et al. Int. Pat. Apl. Pub. No. 
WO 95/15147 ("Ennen"). 

Applicant does not concede that any of the statements by the Patent Office regarding 
what is taught by Ennen are accurate, as the Patent Office has not provided an English translation 
of Ennen. However, Applicant submits herewith a translation of Ermen (Appendix 4), obtained 
in response to the rejections by the Patent Office in its latest Office Action. Ennen appears to be 
directed to the treatment of neurosensory phenomena. Ennen does not teach or suggest a hostile 
biophysical environment comprising an ionic salt mixture comprising choline chloride, sodium 
chloride, and magnesium chloride, as is recited in claim 33, nor does Ennen teach or suggest 
packaging selected from the group consisting of a liposome, an emulsion of collagen, and a 
collagen peptide. Thus, it is believed that claims 33 and 61 are not anticipated or rendered 
obvious by Ennen, and it is respectfully requested that the rejection of these claims be 
withdrawn. 

CONCLUSION 

In view of the foregoing remarks, this application should now be in condition for 
allowance. A notice to this effect is respectfully requested. If the Examiner believes, after the 
foregoing remarks, that the application is not in condition for allowance, the Examiner is 
requested to call the undersigned at the telephone number listed below. 

If this response is not considered timely filed and if a request for an extension of time is 
otherwise absent. Applicant hereby requests any necessary extension of time. If there is a fee 
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fal ulcers and hypertension in salt 



^am Haldiya***, Murli L. Mathur*, 
^^achdev^ and Habibulla N. Saiyed^ 

Jif icine Research Centre (ICMR), Jodhpur 342 005. India 
Mfnstitute of Occupational Health (ICMR), Meghan! Nagar, 
ffSd/^JSO 816, India 

^rpcess of salt manufacture, brine rich in salt is 
^|broad pans exposed to direct sunlight. Brine 
Tslkeep their feet in the brine while working in 
^ins^ They frequently have traumatic ulcers on 
P^and hands. The aim of the present communica- 
la^fto find out if prevalence of hypertension and 
Rfessure in brine workers was affected by the 
^l^f ulcers on their limbs. In our cross-sectional 
^pibnal study, all workers (n = 218) were clini- 
Kltnined. Blood pressure was measured. Systolic 
^stolic blood pressure of brine workers having 
fjbn their limbs was compared with those not 
ggljthe ulcers. In brine workers having ulcers, 
^pnce of hypertension was 15.6%, while it was 
Hn those not having ulcer. Prevalence of hyper- 
^r§was significantly higher in the group of brine 
l^rls with ulcers (P = 0.036) and this may be due to 
yptlbii of salt through the damaged skin. 

Manufacturing is an important industry in Rajasthan, 
^fig more than 20,000 workers*. In the process of 
n^nufacture, highly concentrated brine (groundwater 
ffi^salt, used for salt manufacture) is filled in broad 
g^de on the surface of land, exposed to direct sun- 
iWater evaporates from this brine and sodium chlo- 
fe#}ystallized at the bottom of the pan. In these pans, 
^vorkers sweep crystals of salt using broad wooden 
^knd make heaps of salt near the edge of the pan. 
Efeet and lower portion of their legs remain submer- 
ffibrine for 4 to 8 h per day. Often their upper limbs 
^0 wet with brine or are covered with fine salt parti- 
phey frequently sustain minor occupational injuries 
|ir feet and hands resulting in the formation of ul- 
|yhich heal slowly. There is a possibility of absorp- 
significant amount of salt through such ulcers, 
^in turn may affect their blood pressure. High salt 
|nt in diet is known to increase blood pressure as 
Mf prevalence of hypertension in salt-sensitive sub- 
faThe aim of the present communication was to find 
^prevalence of hypertension and blood pressure of 
^^orkers was affected by the presence of ulcers on 
Jlimbs, 

iSK was a cross-sectional observational study. In Ra- 
bn, occupational health check-up camps were held 
^ree salt-manufacturing sites, under a project on 

Mprrespondence. (e-mail: haldiyakr@rediffmail.com) 
S|eNT science, vol. 87, NO. 8, 25 OCTOBER 2004 



'Prevention and Control of Occupational Health Hazards 
Among the Salt Workers*, sponsored by the Ministry of 
Health, Government of India. This project was approved 
by the Scientific Advisory Committee of National Institute 
of Occupational Health, Ahmedabad, India. The camps 
were organized at Sambhar, Nawa and Phalodi in colla- 
boration with owners of salt-manufacturing units and De- 
partment of Salt, Government of India. Each camp was of 
five days duration. All the workers from nearby salt- 
manufacturing units were invited for their free health exa- 
mination. The workers who were absent on the dates of 
the health camp were excluded from the study. All work- 
ers whose nature of job involved working in brine pans, 
were classified as brine workers. Brine workers who were 
also engaged in other processes related to salt like loading, 
weighing, transporting, milling or packing of salt were 
excluded from the analysis. 

Eight hundred ninety-one salt workers attended the 
camps, out of which 218 who fulfilled the above inclusion 
criteria were taken as brine w'orkers for the present analy- 
sis. Their age in completed years, sex, detailed occupa- 
tional history, including exact nature of job and duration of 
working in salt industry were recorded in performa spe- 
cially designed for occupational health examination. 

After obtaining their informed consent, clinical exami- 
nation including detailed dermatological examination was 
carried out by one of the authors, who did not measure 
blood pressure. Blood pressure of all brine workers was 
measured in supine position after rest of five minutes. 
Blood pressure was measured in the right arm using digi- 
tal blood pressure equipment (Omron T-4). The cuff size 
was 25 cm x 13 cm. Three readings were taken by. trained 
field investigators under the supervision of the first author. 
The first two readings were to familiarize the subjects 
with the process and the third reading was taken for 
analysis. The field investigators were trained to measure 
blood pressure for fifteen days by the authors. Body 
weight and height were measured by another field inves- 



Table 


I. Characteristics 


of study subjects 






Brine workers 


Brine workers not 






having ulcer/s on 


having any ulcer 




Characteristics 


limb/s (;i = 32) 


on limbs (/i = 186) 


P-value 


Age (years) 


29.72 ± 9.78 


31.0418.87 


0.10* 


Gender M/F (%) 


93.7/6.3 


84.4/15.6 


0.16** 


Literacy (%) 


37.5 


33,3 


0.65** 


Income (Rs per 


23740 ± 10076 


23456114615 


0.10* 


annum) 








Smokers (%) 


56.3 


46.8 


0.32** 


Alcohol users {%) 


9.4 


8.1 


0.80** 


BMI kg/m^ 


18,6912.06 


18.81 12.14 


0.10* 


Duration of salt work 


10.63 ±9.17 


10.91 ±7.03 


0.10* 


(years) 








Vegetarians (%) 


65.6 


64.5 


0.90** 



*Studenfs / test was applied. 
**Chi-square test was applied. 
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Table 2. Prevalence of hypertension (HT) according to age 



Brine workers having ulcers Brine workers not having ulcers 



Age in years 


No. examined 


Prevalence of HT (%) 


No. examined 


Prevalence of HT (%) 


15-19 


4 


0.0 


17 


0.0 


20-29 


13 


15.4 


62 


4.8 


30-39 


5 


20.0 


68 


2.9 


40-49 


9 


22.2 


32 


3.1 


50+ 


1 


0.0 


7 


28.6 


Total 


32 


15.6* 


186 


4.3* 



*r = 4.39, P = 0.036. 

tigator trained for the purpose. Height was measured in 
centimetres, using anthropometric rod, while the subject 
stood erect on a flat platform. Body mass index (BMI) 
was calculated as [(weight in kg/height in metres)^]. Sys- 
tolic and diastolic blood pressure of brine workers having 
ulcers on their limbs was compared with those not having 
ulcers. Hypertension was defined as systolic blood pres- 
sure more than 139 and/or diastolic blood pressure 90 or 
above. Student's t test and chi-square test were used to 
study statistical significance of the differences. 

Out of 218 brine workers, 32 (14.7^o) had ulcer/s on 
their hands, feet or lower portion of their legs. The chara- 
cteristics of the study subjects are depicted in Table 1 and 
were comparable in brine workers with ulcers and those 
without ulcers. Both groups did not have significant differ- 
ence in prevalence of smoking, alcohol use, age, literacy, 
income, duration of working in salt industry and BMI. 

Prevalence of hypertension was significantly higher in 
those with ulcer (15.6%) than in those not having ulcer/s 
(4.3%) (;ir^ = 4.39, P = 0.036; odds ratio = 4.12). All 
hypertensive workers with ulcers had systolic hyper- 
tension, but diastolic blood pressure in all of them was 
less than 90 nun of mercury. 

Mean systolic blood pressure in thosq having ulcers 
was 121.84 ± 12.42 mm of mercury and in those not hav- 
ing ulcers was 117.73 ± 12.13 mm of mercury. The dif- 
ference was statistically not significant (/= 1.77, df =216, 
P = 0.08). Mean diastolic blood pressure in those having 
ulcers was 67.44 ± 10.65 mm of mercury and in those not 
having ulcers was 69.11 ± 8,82 mm of mercury; the dif- 
ference was not statistically significant (/ = 0.96, df = 
216, P> 0.10). 

In the present study, overall prevalence of hypertension 
among brine workers was 6.0%, which was higher in those 
having ulcers (15.6%; Table 2). Higher prevalence of hy- 
pertension was consistently observed in all age groups in 
workers having ulcers, while among workers not having 
ulcers, prevalence of hypertension was high only in those 
above 50 years of age. However, hypertension has not 
been found to be an occupational health problem of rock 
salt-mine^ and harbour workers'*. Overall prevalence of 
hypertension in rural population of northwest India^'^ has 
been found to be 3.4-7.2%. 



In India, salt is mainly manufactured from :tW 
shore and subsoil brine. Rajasthan is the second;]!;, 
producer of salt from subsoil brine in India. Salt wo 
are low-paid and do not use gum boots or othelJ 
measures. As a result, their skin is often damaged,^ 
crystals and the spade used for sweeping thern|^S 
workers work in pans containing concentrated brijif' 
their feet and lower legs remain submerged in brin^ 
present study depicts higher prevalence of hypelten 
among brine workers having ulcers on their iini| 
probable explanation of this observation may be'Hlp. 
sibility of absorption of salt through the skin. Wtie 
concentration in fluid surrounding the skin is higWK 
tic pressure and ion-exchange make a significaipt^ 
bution to flux enhancement by diffusion ppte|[ 
Presence of ulcers on the skin is likely to increase!? 
sorption as percutaneous absorption of vario.^ 
stances, including salt is known to be highei^J 
damaged skin*°~*^. This might cause elevation-gj| 
pressure in brine workers having damaged skin0' 
this hypothesis needs to be confirmed. Zhou et^^a^ 
ied the relationship of dietary patterns with Woq(|^ 
in ten population groups and found mean daily^jmj^ 
sodium showed significant positive association):^^ 
systolic and diastolic blood pressure in all groug|§& 
the fishermen. This observation supports our h^ 
as percutaneous absorption of varying quantili|^ 
among fishermen may be the cause for lack of a|^j 
between dietary salt intake and blood pressuix^n^ 
The limitation of the study is that the workers.^.e 
followed and it could not be known whether hy|||t 
disappears after healing of ulcers or disconti|u|tl 
work for few days. Furthermore, biochemical |n|g 
tion of serum and urinary electrolytes could Wjf^^ 
due to lack of facility. Frequent damage of skin^ 
sequent healing may be an important risk factor| 
velopment of . hypertension in salt workers. FurtKerj 
are needed to elucidate our findings. ^ 



I. Haldiya, K. R., Mathur, M. L., Sachdeva, R., Beni||! 
Singh, M. B.. Yadav. S. P. and Lakshminarayana, /J^;J^ 
pattern of desert population engaged in salt manufacjure 
sthan. y. Indian Med, Asxoc, 1995. 93, 95-97. 



Best Available Copy 



06. F. et «/., Dietary patterns in 10 groups and the relation- 
"wfth blood pressure. Collaborative Study Group for Cardio- 
Diseases and Their Risk Factors. Chin, Med. /, 1989, 

«?uliivan, J. J. and Parker, G. D., Investigation of the blood 
"^re levels of workers occupationally exposed to salt. Occup. 

tfoiija, A. and Kontosic, I., Arterial hypertension in correla- 
vi/jth age and body mass index in some occupational groups in 
h|i?li)0ur of Rijeka, Croatia, hid. Health, 1998, 36, 312-317. 
ISotra, P.. Kumari, S., Kumar, R., Jain, S. and Sharma, B. K., 
^filinct and determinants of hypertension in an un-industria- 
.Jfrural population of North India. 7. Hum. Hypertens,, 1999, 
'||7if472. 

i^Ji^N' K. and Kaur, P.. Role of various risk factors in the epide- 
m^$y of hypertension in a rural community of Varanasi district. 
® ^- ^w^'''^ Health, 1 996, 40, 7 1-76. 
mari' P. and Chaudhary, V., Epidemiological study of hyperten- 
giSih^a rural community of western Rajasthan. Indian Heart 7., 
"^3,43-45, 

iglJ^y. N., Kalantri, S. P., Gupta, O. P., Jain. A. P. and Gupta, 
.;^The prevalence of hypertension in rural population around 
ewgram. Assoc. Physicians India, 1993, 41, 422-424. 
'^oncn, J., Murtomaki. L. and Kontturi, K., Effect of diffusion 
gnlial, osmosis and ion-exchange on transdermal. drug delivery: 
3y and experiments. J. Control Release, 1998, 56, 33-39. 
o:dn, -K. C, Wester, R. C. and Maibach, H. I., Diseased skin 
Iwle&in the hairless guinea pig: in vivo percutaneous absorption. 
wiatologica, 1 990. 1 80, 8- 1 2. 
^augh, R. L. and Stewart, R. R, Methods for in vitro percuta- 
^S: absorption studies V: Permeation through damaged skin. /. 
l^.Sci., 1985. 74. 1062-1066. 
^fnirrJ. et aL, Skin penetration of minerals in psoriatics and 
uineaTpigs bathing in hypertonic salt solutions. Pharmacol, Res. 
'"Iwm/i., 1985. 17. 501-512. 

^P:EDGEMENTS. Ministry of Health and Family Welfare, 
^^India had fmanced the project on salt workers working in de- 
!;^P Gujarat and western Rajasthan. We thank the Salt Com- 
er. \Government of India and his team for cooperation. 

a|12 March 2004; revised accepted 1 4 June 2004 



Iflir Ghaggar, Saraswati? 
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identity of the river along which the famous 
a^gan Civilization developed and the causes of the 
||of this culture are topics of considerable de- 
^^any of the Harappan sites are located along 
Ifihenieral Ghaggar river within the Thar Desert 
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in the northwestern India and adjacent I^akistan. The 
Ghaggar was also thought to be the mythical river 
Saraswati, which was described as glacial-fed river. Sr 
and Nd isotopic composition of the Ghaggar alluvium 
as well as Thar Desert sediments suggests a Sub-Hima- 
layan sediment source, with no contribution from the 
glaciated regions. The development of extensive 
Harappan Civilization all along the Ghaggar suggests 
a catchment with high monsoon rainfall. It is likely 
that with the changes in the monsoon scenario after 
3500 BC could have gradually dried up the Ghaggar 
river and resulted in the migration and/or extinction 
of the Harappan Civilization on this river. 

The largest and the oldest urban civilization of the world 
was the Indus Valley (Harappan) Civilization of north- 
west India and Pakistani Almost two-thirds of nearly 
1500 archaeological sites of this civilization occur on the 
dried banks of the Ghaggar river^ (Figure 1). The River 
Ghaggar originating in the Sub-Himalayas flows through 
the northern part of the Thar today as an ephemeral river 
mainly during the SW monsoon season and disappears in 
the desert. However, the river seems to have played a key 
role in the development of the Harappans^'^. The Ghaggar 
river has been identified with the mighty glacial-fed river 
Saraswati"*^, which is described in the oldest religious 
document written in Sanskrit, the Rig-Veda (1500 BC)^'^. 
Based on geomorphological studies and identification of 
clasts in the river channels of outer Himalayas, it has 
been suggested that the palaeo-Ghaggar (alias Saraswati) 
had its catchment in the glaciated Higher Himalayas^'*^. 
Another prevalent hypothesis is that the ancestral channels 
of Yamuna and Satluj once fed the Saraswati^'*'. The an- 
tecedent Yamuna and Satluj rivers originate from the gla- 
ciated Higher and Tibetan Himalayas respectively, and 
limit the expansion of the Thar Desert in the east and 
north. If water availability is the key climate determinant 
for life^^ and the region was already experiencing ari- 
dity'^, the palaeo-Ghaggar must have been perennial for 
the Harappans to flourish. To identify the source (glaciated 
or non-glaciated terrains) and the nature (perennial or 
ephemeral) of the palaeo-Ghaggar and, therefore, to un- 
derstand the likely cause of the social collapse of the 
Harappans, we have studied Sr and Nd isotopic character- 
istics of the sediments deposited by desert-forming pro- 
cesses and by the River Ghaggar in the Thar Desert 
region of northwestern India and compared these with the 
values of other Himalayan rivers in the region and with 
those of various Himalayan lithotectonic units. 

Recycled sediments formed due to cannibalistic pro- 
cesses of erosion and subsequent deposition pose diffi- 
culties for the identification of their sources. However, in 
the Himalayan orogen, Sr and Nd isotopes have been suc- 
cessfully used to identify different lithostratigraphic 
zones and to locate the source areas of sediments of the 
Tertiary foreland basins'^'*^, Bengal fan*^ and those of the 
present-day Himalayan rivers*^. Identification of sediment 
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Abstract 

The permeability characteristics of an organotypic epidermal culture model derived from rat epidermal keratinocytes, ROC, and isolated 
human cadaver epidermis, HEM, were compared. Rat epidermal keratinocyte (REK) ceU line was grown for. 3 weeks on colleen gel m 
the absence of feeder cells in culture inserts at an air-Hquid interface. Transdermal permeabilities of 18 compounds rangmg from 92 to 
504 in molecular weight and from -4.3 to 3.9 in log. of octanol-water partition coefficient, charged or uncharged, were measured m the 
culture model and isolated human epidermis. The REK organotypic culture model (ROC) provided a close ^timate of human epidermal 
oermeabiUties over the whole range of the solutes used with on the average of 2.fold higher permeability coefBcrents (range 0.3-5.2) than 
fose obtained from isolated'hM^ The easily inaintained and ireprodudble ROC iu^^^^^ be useM m screemng 

transepidermal drug peniieabilities together with possessing potential for research on dermal formulauons, irritation, toxicity and gene 

therapy. . - ; 

© 2003 Elsevier B.V. All rights reserved. 
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1. Introduction 

Aciequate skin permeation is a prerequisite for success- 
ful (trans)dernial drug therapy. Isolated human or animal 
skin is c^^iiimionly used in drug permeability assessment: 
The problems involved with isolated skin samples include , 
"TEeir* variable quality and lack of viability. 

It is widely acknowledged that organotypic epidermal 
culture models may be used to replace the human cadaver 
and animal skin that are difficult to obtain and be stan- 
dardized for routine permeability testing of drugs. Several 
skin cell culture models have been presented in the 
hterature.* Typically, they are either derived from primary 
or early passage ceUs,.e.g. (Fartasch and Ponec. 1994; 
Kennedy et al., 1996; Stark et al.. 1999) or from continu- 
ous ceH lines (Boelsma et al., 1999; Schoop et al., 1999). 
These cultures are generally supported by feeder cells 
(AsbiU et al., 2000; Boyce and Williams, 1993; PouHot et 
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al., 1999; Schoop et al., 1999) or de-epidermized dermis 
(Cumpstone et al., 1989; Fartasch and Ponec, 1994; 
Kennedy et al., 1996; Mak et.al., 1991). . 

Keratinocytes maintained at the air-liquid interface have 
been shown to stratify, express several differentiation 
markers, and form an apparentiy normal stratum comeum 
(Fartasch and Ponec, 1994; Mak et al., 1991; Parenteau et 
al., 1991; Schoop et al., 1999; Slivka et al., 1993). 
However, it has been difficult to produce an organotypic 
culture model with a permeability barrier comparable to 
that of intact skin (AsbiU et al., 2000; Cumpstone et al., 
1989; Mak et al.. 1991; Nolte et al., 1993; Ponec et al., 
1990; Regnier et al.. 1993; Roy et al., 1993; Schmook et 
al.. 2001). Typically, the permeabilities .of the cultured 
skin alternatives, derived either from human or animal 
skin, have been several orders of magnitude higher tiian 
that of intact human skin in vitro. In. addition, so far no 
extensive evaluation of drug permeabilities with a large set 
of compounds with a variety of physicochemical charac- 
teristics in any cultured epidennal models has been done. 
This is obviously a serious defect, in the use of cultured 
epidermis if they are to be used to predict human skin 
pejmeabilify, and may lead to erroneous conclusions. 
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The purpose of the present study was to evaluate a novel 
permeability model, which', at the same time, would be 
relatively simple and easily maintained. Recently, an 
organotypic epidermal culture model utilizing a continuous 
keratinocyte cell line derived from newbom rat skin (REK) 
was shown to differentiate into epidermis with normal 
stratum comeum ultrastructure and functional permeability 
barrier as measured by trajnsepidermal water loss (TEWL) 
and by permeability of corticosterone (Pasonen-Seppanen 
et al., 2001a,b). Namely, after 3 weeks in culture, the 
transepidermal water loss values and the permeability 
coefficient for corticosterone, a lipophilic solute, were 
close to those measured in intact human cadaver epi- 
dermis. The present work further explores the permeability 
barrier properties of the REK organotypic culture (ROC) 
model. It is shown for the first time that the permeability 
coefficients of ROC for a large number of solutes with 
different Upophilicities, molecular weights, and charges are 
close to those measured in intact human cadaver skin in 
vitro. 

2. Materials and methods 

2.1. Materials 

Cell culture media. Minimal essential medium (MEM, 
without L-glutamine) and Dulbecco's MEM (DMEM, with 
4500 mg/1 glucose) were obtained from Gibco-BRL (Life. 
Technologies, Paisley, Scofland), as well as Earle's Bal- 
anced salt solution (EBSS, lOX), and 7.5% sodium 
bicarbonate solution. Fetal boyine serum (FBS) was ob- 
tained from HyClone (Logan, UT, USA). L-Glutamine 
(200 mM), penicillin-streptomycin solution (containing 
10 000 U penicillin and 10 mg/ml streptomycin in 0.9% 
sodium chloride), trypsin-EDTA solution (containing 500 
U/ml porcine trypsin and 180 ixg/ml EDTA in Dulbecco's 
phosphate-buffered saline), and L-ascorbic acid were from 
Sigma (St. Louis, MO, USA). Transwell tissue culture 
inserts (24 nam diameter, 3.0 |xm* pore, size) were from 
Costar (Cambridge, MA, USA). Aqueous acetic acid and 
sodium chloride were of standard laboratory grade. 

For permeation studies, phosphate-buffered saline (pH 
7.4, prepared from PBS tablets, Sigma) was used. The 
radiolabeled permeants (fructose-l-^H(AO)-sucrose (12.3 
Ci/namol), i>-(l-^H(A^))-manmtol (26.3 Ci/nunol), 3-0- 
(methyl-^H)-methyl-D-glucose (81.5 Ci/mmol), (2-^H)- 
glycerol (16.0 Ci/nmiol), ■D-(l,2,6,7-^H(iV))-aldosterone 
(70.0 Ci/mmol), (l,2,6,7-^H(i\r))-corticosterone (70.0 Ci/ 
mmol). (la,2a-^H(iV))-testosterone (51.0 Ci/mumol), (6,7- 
^H(A^))-estradiol (44.0 Ci/nmiol), and (7-^^C)-saKcylic 
acid (55.5 Ci/mmol) were obtained firom^ NEN .Life 
Science Products (Boston, MA, USA), and ^H(G)-raffinose 
(20.0 Ci/mmol) from American Radiolabeled Chemicals 
(St. Louis, MO, USA), all at higher than 97% purity 
according to the manufacturers and used as such. The 
P-blocking agents atenolol, sotalol . hydrochloride,, nadolol, 
pindolol, metoprolol tartrate, alprenolol hydrochloride; and 



propranolol hydrochloride were obtained from Sigma, 
except for timolol maleate, which was donated by Merck, 
Sharp & Dohme Research Laboratory (Rahway, NJ, USA), 
all at higher than 98% purity and used as such. The 
physicochemical properties of the solutes used are pre- 
sented in Table 1. 

2.2, Keratinocyte culture 

The rat epidermal keratinocyte cell line (REK) was 
derived from neonatal rat epidermal keratinqcytes (a 
generous gift from Dr. Donald K. MacCaUum, Ann Arbor, 
MI, USA) ^originally isolated by Baden and Kubilus 
(Baden and Kubilus, 1983). Stock cultures were grown in 
MEM vidth 10% fetal bovine serum, 4 mM L-glutamine, 50 
|jLg/ml streptomycin sulfate and 50 U/ml penicillin at_ 
37 °C in humddified 95% air/5% COj.. The REKs' were 
subcultured twice a week by iucubating them for /Approxi- 
mately 5 min at 37 °C in trypsin-EDTA solution. The 
REK cell line was tested to exclude the possibility of 
Mycoplasma contamination. r ' 

2.3, Collagen substrate 

Approximately 1 g of collagen strands from three to four ' 
rat tails was removed, washed with hQ% ethanol and 
allowed to dry under UV light. Dry collagen strands were 
weighed and sterile aqueous acetic acid (0.1%, v/v) was 
added to give a final concentration of 3 mg/ml. The 
mixture was stirred for 48 h at 4 °C and then centrifaged 
for 2 h at 2000Xg. The clear supernatant was stored at 
4 °C. To prepare collagen supports, the cqllagen was mixed 
with EBSS, 7.5% sodimn bicarbonate, and 1 .M sodium 
hydroxide solution, at a volume ratio of 80:10:3:2, respec- 
tively, in an ice bath. To cover the entire growth surface, 
1 ml of the mixture was pipetted onto individual tissue 
culture inserts and incubated overnight at. 37°C in a 
hmnidified atmosphere to polymerize collagen, and rinsed 
briefly with the culture medium before use. 

2.4, REK organotypic culture (ROC) 

Recenfly confluent cultures of REICs with no mor- 
phological evidence of stratification were trypsinized, 
diluted witii DMEM with 10% FBS, 4 mM L-glutamine, 
50 |JL,g/ml streptomycin sulfate and 50 U/ml penicillin, . 
and 2 nal of the cell suspension (200 000 cells /ml) was 
applied onto the collagen substrates on the transwell 
culture inserts. The subcultiVated REKs were grown for 3 
days with culture medium present both in the well beneath 
the insert as well as on the surface of the cells. The culture 
was subsequentiy air-lifted (i.e., medium removed from the 
surface of the cells) and the level of die ihedium beneath 
the insert adjusted to the level of REKs present on the 
collagen gels (approximately 1.5 ml). The culture medium 
was supplemented with 40 |Jig/ml ascorbic acid from the 
day after the- culture was configured at the. air-liquid 
interface. The cultures were grown for 3 weeks, with- the 
medium changed every 2 days for the first week and daily 
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Table 1 

The physicochemical and penrieation properties of the solutes used in the present study 



Solute 



MW Cadaver epidemiis 
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0.1910.09 
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5 


0.8 


-1.06" 
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0.2210.09 
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0.2210,11 
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5' 


1.0 


-0.33" 
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0.5510.23 
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0.6410.37 
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0.8310.33 
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0-6' 


10 


2.0 • 


1.94^ 
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Raffinose 

Sucrose 

Mannitol 

Methylglucose ; 

Salicylic .acid 

Glycerol 

Atenolol 

Sotalol 

Nadolol 

Pindolol 

Metoprolol 

Tunolol 

Mprenolol 

Aldosterone 



Corticosterone 

Testosterone 

p-EstradioI 



' Log of octanoi-water partition coefficient • ; 

" Molecular weight. 
' From Yoneto et al. (1995). 
From Johnson et al. (1997). 
"From Hansch and Xeo (1979). 

'From Hanschetal. (1995). = • 

From deGier (1993). . . . ■ 

" Log of octanol/.pH 7.4 buffer partition coefficienL Calculated according to Schbenwald. and Huang (1983); Wang et:al, (19^1). r ; , 
^ From Johnson et al. (1996). 

'From Yallcowsky etal, (1983). . . . . ' 

' Permeability coeffipient (mean! S.p.). . . 

" Lag time of permeaiion. u \ 

" Number of individual samples (for cadaver sldn, samples were from several donors; for culture, samples were of at least two passage numbers). 
* Permeability data obtained from the literature as referenced in Ref . column. 

P The permeahility coefficient of cultured epidermis divided by the permeabiUty coefficient of cadaver epidermis. 



thereafter. Histology and ultrastructure of the cultures were 
evaluated periodically in a manner similar described in 
previous reports from our group (Pasonen-Seppanen et al., 
2001a,b)y 

''IS" Permeation studies 

For the permeation experiments, 3-week-old REK or- 
ganotypic cultures (ROC) or human epidermal membrane 
(HEM) was used. Cadaver abdominal skin was obtained 
from the Kuopio University Hospital (Kuopio, Finland) 
with the permission from The National Board of 
Medicolegal Affairs. The excised skin was heated in 
distilled water at 60 °C for 2 min before the epidermis was 
separated fomi the underlying dermis, dried, and frozen for 
later use. Prior to mounting the pieces of skin into 
diffusion chambers, the frozen skin samples were allowed 
to thaw at room temperature, and hydrated briefly in PBS 
prior to mounting the skin in the diffusion cells. 

The permeation studies for both HEM and ROC were 
carried out using two-chamber diffusion cells (Side-Bi- 
Side, Crown Glass Company, Somerville, NJ, USA) with a 
3-ml volume in each compartment and an average effective 



divisional area of 0.64 cm^ Briefly, the membrane was 
clamped between the cell-halves with the stratum comeum 
side facing the donor chamber. HEM was hydrated Qvqr- 
night in PBS at 3TC prior to the experiment Continuous 
magnetic stirring was maiutained in both chambers. At the 
beginning of an experiment, the receiver chamber was 
replaced with fresh buffer. The donor, chamber was 
replaced with a solution, cpntaining 1 mM of each of the 
p-block6rs studied in PBS (pH 7.4). Altematively, the 
donor was spiked with a radiolabeled test substance 
(20 000-40 000 dprn/lO fxl). The concentrations used 
were far below the solubility limits of the solutes. Aliquots 
were withdrawn at predeterinined time intervals (typically, 
10 |xl from the donor and 200 fxl from , the receiver 
chamber). The same volume of fresh buffer, was added 
back to the receiver to maintain a cpnstant volume. The 
samples were analyzed for the (3-blocking agents by HPLC 
(see below). The radiolabeled samples were mixed with .3 
ml of scintillation cocktail (UltimaGold, Packard, Bio- 
science^ Groningen, The Netherlands) and analyzed by a 
liquid scintillation counter . (WinSpectral 1414, Wallac, 
Fiol^nd). 

. Permeability experiments were run long enough to 
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ensure .that the steady-state phase was clearly reached and 
maintained (typically, total time for data collection was at 
least three times the lag time). The permeability coeffi- 
cients (P, cm s~^) for the probe permeants were calculated 
at steady state under sink conditions according to Eq. (1): 



P = 



1 dQ 



(1) 



where A is the diffasional area of the diffusion cell (cm ), 
is the concentration in the donor chamber (dpm ml ~ 
or jxg ml"^), and dQ/dt (dpm s~^ cm"^ or |xg s"^ cm"^) 
is the slope of the linear region of the cimiulative amount 
of permeant .in receiver chamber ((2, dpm or fig) versus 
time (?, s) plot with the values being in the range of 
0.941-1.000, except for the slow penetrating permeants 
such as atenolol and sucrose where r^= 0.799-1.000. 

The integrity of the skin samples, and presence of air 
bubbles were monitored visually throughout the study. The 
data from chambers showing a decrease in medium 
volume, a significant decrease (i.e., more than 5%) in 
donor concentration, or an abrupt increase in the steady- 
state slope, were discarded; 

The experiinental data obtained with ROC and HEM for 
the 18 compounds studied was fitted with Eq. (2): 



log P (em s"') = a • log^fiT,,, + ^ • MW + y 



(2) 



where P is skin permeability (cm s"^), 
water partition coefficient, MW is molecular 'weight, and 
a, fi, and y are constants. The multiple regression analysis 
was- iperformed lising the SigmaPlot 4.01 software (SPSS 
Science, Chicago, IL, USA).. • 

a ■ 

2.6: :H]^ assay 

A Beckman Gold HPLC system (Berkeley. CA, USA) 
with a Kromasil 100 C8 column (5 ixm, 4.6 mmX15 cm, 
Higgins Analytical, Mountain View, CA, US A), was used 
to determine the ^-blocker concentrations in the samples at 
30 T using gradient HPLC ais described previously (Ranta 
et al., 2002). Calibration curves were constructed for each 
experimental run. The peak area was linearly related to 
concentration for samples containing 2-1400 rig of solutes 
injected onto the column. The correlation coefficients for 
the calibration curves were more than 0.98. The peak areas 
of the samples were converted to concentrations using the 
calibration curve. 



K^^^ is octanol- 



3. Results and discuission 

As recently reported, the 3-week-old ROCs sup- 
plemented with vitamin C in the present study were shown 
to duplicate most of the structural features of epiderinis in 
vivo (Pasonen-Seppanen et al., 2001a,b), . Histology and 
ultrastructure of the cultures in the present study were 
checked periodically. The structural features remained 
iinchanged throughout the study (data not shown). 



The penneabUity characteristics of the ROC remained 
fairly constant for almost 60 passages (about 30 weeks) in 
our laboratory, as studied by mannitol (Fig. la), a hydro- 
philic solute, and corticosterone (Fig. lb), a lipophilic 
model solute. Permeability coefficients for these model 
solutes were checked periodically and used as a means to ^ 
check the validity of the culture in the course of various 
experiments. 

It is notei that the present experimental conditions 
included an overnight prehydration of the HEM prior to 
the study to ensure that the water content of the tissue 
would be comparable to that of the cultured epidennis (i.e., 
which is maintained in constant contact with the culture 
medium, and therefore, is- considered to be fully hydrated). 
The long prehydration is not generally included in permea- 
tion experiments with skin. Ho\yever, in generai;'^tBe 
results from permeability experiments with HEM were 
found to be in line with the values reported in the literature 
(again, which mostiy includes results from studies that, are 
conducted without prehydration) (Table 1). Therefore, 
prehydration of the HEM is not considered to play a major 
role in altering skin permeability. Only the permeabihty 
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Fig. 1. Stability of the epidermal barrier properties of ROC as indicated 
by the penneabUity for (a) mannitol and <b) corticosterone in cultures 
raised from consecutive passages of the cell line. 
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coefficient for sucrose obtained in the present study was 17 
times higher than that reported earlier (Peck et.al., 1998), 
the reason for which remains unclear. A possible explana- 
tion may be found firom the non-Gaussian distribution of 
human skin permeabilities (Williams et al., 1992). 

As illustrated in Fig. 2,. the permeability of ROC for the 
various solutes mimics closely that of HEM, There is a 
general increase in permeability with an increase in 
octanol-water partition coefficient in both model mem- 
branes. However, the ROC permeability seems to be more 
sensitive to changes in lipophilicity of the model solutes. 
With the exception of mannitol, it is shown that below 
log K.of approximately -2.0, the permeability of ROC for 
the model solutes underestimates HEM permeability, 
whereas at above log JS* of O.Q, the permeability of ROC is 
usually somewhat higher than that of HEM (with again one 
outlier, aldosterone, for which the permeabilities of ROC 
and HEM are ahnost equal). 

To further explore the selectivity of ROC membrane and- 
to make a coniparison to that oif HEM, the experimental 
data obtained in the present study, was fitted with Eq. (2). 
Eq. (2) describes skin permeability in terms of the 
octanol-water partition coefficient (log K^^^ and molecu- 
lar weight (MW) (a, A and y being constants) and it has 
been successfully used to fit skin permeability data firom a 
large- variety of compounds (Moss and Cronin, 2002; Potts 
and Guyv 1992), The resulting quantitative structure-per- 
meabiUty relationships (QSPR) for the compounds studied 
in ROC and HEM are presented in Table 2 and in Fig. 3. A 
multiple regression' on log jST^ct ^nd MW. according to Eq. 
(2) provided a fit with r^ = 0.48 for HEM and r^.=0.51 for 
ROC. The values of the fitted parameters for HEM were 
found to be somewhat different ' than those found with, 
human skin in the literature (Moss and Cronin, 2002; Potts 
-and Guy, 1992), which may also indicate that the small 
. number of compounds does not allow for accurate model- 
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Fig. 2. The permeability coefficients (as meaii±S.D.; n = 3-21) of human 
cadaver epidermal membrane (HEM) and the cultured skin model (ROC) 
for the various. solutes used in the smdy as a function of their octanol- 
water partition coefficients (log K). 



ing. of the data. In addition, a large number of the 
compounds used in the present study are not included in 
the models presented previously in the literature. 

In general, not a very clear difference was found 
between the obtained values of the fitted parameters for 
ROC and HEM (Table 2). The a value reflects differences 
between the two partitioning domains, i.e., octanol (a -I) 
and skin. The values lower than unity obtained for both of 
the membranes suggest that the partitioning domains of 
ROC and HEM are naore polar than octanol, which 
parallels the findings firom earlier studies (Anderson et al., 
1988). The value of a for ROC (0.25±0,07) was found to 
be somewhat— although not significantly — ^larger than that 
for HEM (0.17±0.05) (Table 2). Consequentiy, ROC may 
be regarded as a slightly less polar partitioning domain 
than HEM. In general, selectivity in solute partitioning 
increases with decreasing polarity of the partitioning 
domain (Anderson et al., 1988), and, consequenay, the 
differences in ROC permeability are greater as the solute 
hpophilicity is varied, which is in. line with the present 
experimental findings (Fig. 2). The reciprocal vaiiie of j3 
may be approximated as a measure of the average firee- 
volume available for diffusion (Guy and Potts, 1992). 
ROC (|iS|=0.0023±0.0013) and HEM (\^\ = 
0.0028±0.0018), in this context, may be regarded equal. 
Sinularly, the absolute value^ pf y (|rl='5.8) for ROC is 
very close to that obtained for human skin (ly| = 6.2), most 
likely reflecting a similar mechanism of transport through 
the tortuous lipid domains of human skin and ROC (Guy 
and Potts, 1992). 

The results suggest possible stracturai or compositional 
differences in the barrier domains of the membranes. 
PrevioTXS work on organotypic keratinocyte cultures shows 
that high permeability exhibited by the cultures probably 
reflects abnormalities in the stratum cdmeum lipid profile 
(Ponec,. 1991; Vicanova et ^il., 1996) as well as the 
intercellular organization of lamellar, body contents (Boyce 
and WilHams, 1993; Fartasch and Ponec, 1994), the main 
determinants of the permeability barrier (Elias and Menon, 
1991).. It is noted that the ROC model is devoid of 
appendages , that may provide an additional route of 
penetration for polar solutes in HEM. In addition, it is of 
interest to explore the lipid composition and structural 
arrangements in the mtercellular lipids of ROC stratum 
comeum, which remain to be elucidated in further studies, 
Despite of the slight differences in selection for lipo- 
philicity among the membranes, there is a good correlation 
(r^= 0.905) between ROC and HEM in terms of overall 
permeability characteristics, as is .evident from Fig. 4: 



Table 2 



ability data using Eq. (2): logP = 



— — O "Oct r 

Membrane 


. a' 


)S*X10' 






HEM 

ROC 

y 


0.17±0.05 
0.25 ±0.07 


-2.3 ±1.3 
.-2.8±1.8- 


-6.20±0.38 18 
. -5.84±0.51 18. 


0.48 
0.51 



* ValueiS.D. derived firom multiple regression analysis. 
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■ y = 2.5062 + 1 .3352x 0.90496 



-4.5 





Fig. 3. A graphical representation of a multiple regression analysis of the 
permeability data obtained in the present study for (a). HEM and (b]rROC 
according to Eq. (2). 



log P (HEM) 

. Fig. 4. The correlation between ROC and HEM permea^^iities for the 
solutes used in the study. 

The present study demonstrates for the first time that an 
organotypic keratinocyte culture may give a. very close 
estimate on HEM permeabilities for solutes over a large 
range of lipophilicities and molecular weights. Although 
the present culture model is derived- from rat keratinocytes 
it shows almost similar permeabilities and a good correla-' 
tion with isolated human epidermis. A model derived from 
human skin cells, of course, would be more desirable, if 
such permeability properties could^be obtained.. However, 
at this point the ROC model seems useful, for drug 
permeability testing. In addition, livcv-organotypic cultures 
also offer potential for research on dermal formulations, 
irritation, toxicity and gene therapy. 
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Importantly, the drug permeabilities in ROC and HEM are 
similar," but not identical, however. Of 18 tested com- 
pounds, 14 showed 3-fold or smaller difference in per- 
meabilities of the- model membranes. The remaining four 
of the 18 compounds exhibited larger than 3-foldj but less 
than 5.2-fpld, difference in penneabilities in the two 
membranes. Taken, together, the performance of the ROC 
" model in predicting human skin permeabilities is reason- 
able. 
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Abstract 

Pick's first law describes diffusive flux as a linear function of the concentra- 
tion gradient; its most popular generalization describes flux as a linear func- 
tion of the chemical potential gradient. This generalization and others have 
been used for modeling, but the relationships between the flux and the gra- 
dients are nonhnear because the coeflSicients generally depend on state vari- 
ables. In this paper, lattice density functional theory equations for diffusive 
flux are recast into linear differential form to explore weil-knovra diffusion 
equations. This shows that the diffusive flux is related to the gradients of state 
variables (e.g., to the chemical potential gradient) through integrating factors, 
much like the heat flux in Clausius's theorem^ is related to the difference in en- 
tropy. Subsequent analyses show why diffusive flux is more hnear with re- 
spect to the fugacity gradient than the chemical potential gradient; and that 
the gradient of another property has linear proportionaUty to the diffusive 
flux. This property can appear as the impingement rate onto vacancies and 
molecules of species whose density gradients can be influenced by diffusion; 
and it behaves similarly to the chemical potential and the fugacity. The re- 
sults agree with the theory of non-equihbrium thermodynamics about the 
conjugate force vectors for the diffusive flux. 

1. Introduction 

The description of diffusive flux in the form of a linear law is used commonly 
to model diffusion. The classic example is Pick's first law, 
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Jn^'DnVCn, (1) 

where n indicates the species,'-/„4S^ the diffusive flux in units of molecules per 
area per' time, Z)„ is the transport diffusivity^andjQiJs the concentration [1- 
3]. Although Eq. (1) has a linear form, the proportionality between the diffu- 
sive flux and the concentration gradient is linear only when the diffusivity is 
constant; under these conditions, the concentration gradient appears to be the 
"driving force" for the process. The diffusivity can be constant for uniform 
systems, such as isothermal and isoviscous phases with a uniform total den- 
sity; it can be predicted by using the Einstein or Stokes-Einstein relation [2, 
4, 5] and by using molecular simulations [6-10]. For diffusion of vapors, the 
partial pressure gradient often is used to describe the diffusive flux [11], 

- -^«VP„, (2) 

where Pn is the permeability, and Pn is the partial pressure of species w. For 
isothermal systems in the ideal-gas limit, the permeability can be written as 
Dl/kT, where the superscript indicates the gas phase and kT is the ther- 
mal energy; and through membranes, the permeabihty can be written as the 
product of solubility and diffusivity. Equation (2) has a linear form, but the 
proportionaUty between the diffusive flux and the partial-pressure gradient is 
linear only when the permeabihty is constant - in such a case, the partial- 
pressure gradient appears to be the driving force for the process. For diffusion 
in nonideal solutions, the fugacity gradient can be used to model diffusion, 

l = -L{Wfn^ (3) 

where L{ is the coefficient associating the diffusive flux and the fugacity gra- 
dient V/„, and where the fugacity is related to the chemical potential by 
— fi^^ = kT \n{fn/Pn)' For isothermal ideal solutions, can be written 
DlCtoi/P^\ where the superscript '7" indicates the Uquid phase, Qot is the 
total concentration, and P^^^ is the saturation vapor pressure of species /t. 
Equation (3) has a linear form, but the proportionality between the diffusive 
flux and the fugacity gradient is linear only when the coefficient is con- 
stant - in such a case, the fugacity gradient appears to be the driving force 
for the process. Recent isothermal studies of multicomponent pervaporation 
indicate that the fugacity is the fundamental variable for describing fluxes 
and separation factors for zeoUte membranes [12], 

The fundamental driving force for diffusion often is assumed to be the 
chemical potential gradient. In the absence of irreversible processes that are 
coupled to diffusion, the diffusive flux can be written as a Unear fxmction of 
the chemical potential gradient. 
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/„ = ->f'^'V//„, (4) 

where ii„ is the chemical potentiaLQC.s|)ecies n, and ^ is the mobility [13]. 
Assuming that the chemical potential graHieht-(i)-dadxejjhffusio^ in general 
and (u) appears as the concentration gradient for ideal systems, while express- 
mg the chemical potential as +/:rhi(y„C„) in Eq. (4), gives the 

Darken-Dehlinger relation [14, 15] for the transport diflFusivity in Eq. (1), 

(5) 

where Z)^ is the self-diffusivity and y„ is the activity coefficient. Accordingly, 
soktions with constant activity coefficients have a mobility Lj^" that equals 
DlCn/kT , which depends on the position along a concentration gradient for 
a constant diffusive flux. The transport diffusivity in Eq. (5) also was derived 
by Onsager and Fuoss [16], and it can be written [17] as 

^"-^"51^' (6) 

where p„ is the number density. For a historical description of Eq. (6), see 
Reyes et al. [18]. That the chemical potential gradient drives diffiision'is a 
paradigm that generally is attributed to Einstein [2, 4]; however, historical 
evidence .suggests that this principle dates back to Maxwell [19]. Nevertheless, 
any description, of the diffusive flux that assumes Eq. (4) departs from Pick's 
first law according to Eqs. (5) and (6), where D° is sometimes approximated 
as being constant [17]; this includes a Maxwell-Stefan formulation using 
D°„=D„, [20]. ■■ 

The strongest theoretical support for Eq. (4) and, consequently, for Eqs. 
(5) and (6), appears in the dissipation-phenomenological equation (DPE) 
approach [21] of linear non-equilibrium thermodynamics (NET) [22], which 
restricts itself to Unear phenomenological laws but admits that "It is very 
well possible that some irreversible processes must be described by non-Unear 
phenomenological laws" [23]. The NET-DPE approach is "thermodynami- 
cally sound", and it can be used to analyze irreversible processes that are 
coupled [24]. NET describes the entropy production due to diffusive transport 
as a summation of net fluxes multiplied by conjugate forces, i.e., JiXf, and 
it describes the net fluxes as Unear functions of these conjugate forces, i.e., 
= J2k^ik^fc> unhke Eqs. (l)-(4), where the coupling of processes is not ob- 
vious. The conjugate forces that multiply the fluxes are sometimes called 
''thermodynamic forces" or "affinities", but this paper refers to them only as 
"conjugate forces". It is known that the phenomenological coefficients Luc 
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generally depend on the state variables, which usually depend on position and 
time; this leads to nonlinearity as well as unequal partitioning of entropy pro- 
duction for equipartitioned forces [25, 26]. To identify the conjugate forces 
for the diffusive flux equatiSnr -one. must describe the total entropy produc- 
tion of which a portion is caused by diffu^fransport; for a multicomponent 
system we have, 



+ •••> (7) 



where cr is the rate of entropy production per volume, Q is the heat flux vec- 
tor, T is the temperature, the conjugate force vectors Xk for diffusive flux ap- 
pear in brackets, the superscript in fil indicates that this chemical potential 
neglects external fields or long-range interactions, and Ft is a force vector 
that can include contributions from an external field or long-range interac- 
tions, and the ellipsis . .") represents entropy productions due to viscous 
forces and chemical reactions. For conservative forces, Fk = -V^^ where (p^ 
is a potential energy. By modifying the heat flux to Q' = Q - Ejt /^Pk in the 
total entropy production expression, as shown by de Groot and Mazur [23], 
the entropy production due to diffusive flow reduces to 

Q' vr 1 ^ 

^= f2 7^2^Jk-V;zjt + ---, (8) 

k=l 

where ^k=/4 + <Pk (de Groot and Mazur write //^ as fi^); in the transition be- 
tween Eqs. (7) and (8), the entropy flux effectively has been transformed from 
Q/r to Q'/r. As mentioned already, this NET-DPE result provides strong 
support for Eq. (4), because the conjugate-force vector in Eq. (8) is the chem- 
ical potential gradient. At the same time, the NET-DPE approach using Eq. 
(7) shows the limitation of Eq. (4) and, consequently, pf Eqs. (5) and (6). 
Only recently, it has been shown in this journal [27] that Onsager's reciprocal 
relations are not always obeyed when the chemical potential gradient is the 
conjugate-force vector. 

Equations (l)-(4) are examples of linear diffusion equations. Equation (4) 
often is considered to be a general law because the chemical potential is 
a fundamental thermodynamic concept, species usually migrate down their 
chemical potential gradients, chemical potential profiles usually are contin- 
uous at a steady state, and the equation is consistent with hnear NET theory. 
Equation (3) has the same attributes as Eq. (4) because, historically, the fu- 
gacity has been defined directly from the chemical potential. On the other 
hand, Eqs. (1) and (2) are considered less general and can be less predictive 
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because their coefficients can be positive or negative, as well as discontinuous; 
e.g., the diflfusivity in Eq. (1) has negative values in the case of species migrat- 
ing to form separate phases that-are purer than the original mixture [28]; and 
the permeability in Eq. (2) passes through a maximuptp^d can be discontin- 
uous when changing the pressure in porous materials [29-32]. Depending on 
the application, one of these equations may be more suitable than the others, 
but Eqs. (l)-(3) are considered to be forms .of Eq. (4). Interestingly, the coef- 
ficient Ll in Eq. (3) is constant more often than ^ in Eq. (4). For example, 
consider isothermal ideal solutions with constant diffusivities exhibiting Fick- 
ian behavior, as well as a binary isothermal system of nonideal spherical mol- 
ecules undergoing equal counterdifiFusion in an external field [33]. 

Linear diffusion equations such as Eqs. (l)-(4) are important because they 
are relatively intuitive and they can minimize computational overhead, e.g., 
when applied to multicomponent systems [34]. iSowever, Hnear equations for 
the diffusive flux have been considered to be approximations to the Maxwell- 
Stefan equations [34], and they usually are written in terms of gradients that 
do not give a universal linear law (i.e., the following gradients are associated 
with coefficients that generally are not constant): the concentration gradient, 
the pressure gradient, the fugacity gradient, the pH gradient [35, 36], the elec- 
trostatic or diffusion potential gradient [37, 38], and the chemical potential 
gradient. This hst can include a multitude of additional references because 
the problem of formulating a general linear law is more than 50 years old; 
therefore, interested readers are referred to the introduction in an excellent 
paper by Merk [39] that mentions relevant work from the early and mid- 
1900s. 

It would be convenient to have a linear diffusion equation in terms of one 
-gradient with a constant coefficient, 

where is the diffusive flow in units of molecules per time, ^ is a constant 
cross-sectional area through which diffusive flow occurs, X unconditionally is 
constant at a value of unity multiplied by the spatial unit used in evaluating 
the gradient, and r„ is the system's property that drives the diffusive flux 
of species n. Additionally, it would be interesting to understand the funda- 
mental nature of the function r„ (e.g., is it a function of state? Why or why 
not?). Lastly, it would be reassuring if Eq. (9) also gave (a) the correct formu- 
lation of the Maxwell-Stefan equations, (b) correct equilibriimi limi ts^ (c) 
Fickian behaviour with correct departures for nonideal systems, (d) reasonable 
non-equilibrium phase behavior, (e) quantitative agreement with molecular 
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dynamics simulations, (f ) a vanishing diffusive flux at the critical point, and 
(g) agreement vv'ith the NET-DPE approach. 

We have developed a lattice 'density . fun^^ theory (LDFT) approach for 
modeUng diffusion [29] as an extension of eqiiJlibFiuni-LDFT [40-48]. The 
LDFT approach gives item (a), as shown in reference [49] with a detailed 
example for a quaternary system. Reference [29] shows that the LDFT 
approach gives items (b) through (d). To reinforce item (d), it can be shown 
[50] that the approach is consistent with the Giacomin-Lebowitz-Marra 
theory [51] of non-equilibrium phase transitions. Lastly, the LDFT approach 
can give item (e) at least for interdiffusing Weeks-Chandler-Anderson fluids 
in an external field [52]. In this paper, we use the LDFT approach to provide 
Eq. (9), to show that Eq. (9) is related to the well-known Eqs. (l)-(4) through 
iiitegrating factors, and to demonstrate items (f ) and (g). 



2. LDFT equations for ihulticomponent systems 

The LDFT equations for mass transfer in multicomponent systems [52] can 
be simmiarized in the following molecular balance at the general lattice site r: 



5Pn 



dt 



(10) 



where -the number density is limited to 0 < < 1, « is the coordination 
number. One-directional molecular flows (in units of molecules per time) 
occur between sites r and s in the lattice reference frame 



Jn 



1 



Pnir) 



1 - 



(11) 



The characteristic time r„ is related to the mean lifetime of «-type molecules 
at site r; the vacancy probability is (1 - X)/?m); the functional defines the 

influence of the reduced potential energy field on the migration of species n, 
from site r to site s, 



^rs ^ fO; Mr) > Er.{s) 

" \-\E„{t) - E„{s)\; En{r)<E„{sy 



(12) 



where 



Mi) = 



£n(J) + (Pnjl) 
kT{i) 



(13) 
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,The term (p^{r) is a particular long-range interaction or an external potential 
energy applied to molecules at site r; and the intermolecular interaction en- 
ergy En{r) is a sum of nearest^neighbor contributions. For a system with only 
first nearest-neighbor interactions, the interaction:.-energy can be written 



X 

^n{r)=Ylenm^Pm{s), (14) 

n s 



where E^m is the interaction energy between nearest-neighbors of species n and 
m. Equation (12) is analogous to the Metropolis Monte-Carlo algorithm, in 
which one molecxile at site r is moved to a vacant site s with a probability 
equal to e^". However, this work is not based on Monte-Carlo dynamics: 
molecular motion occurs throughout the system; and the terms e^", and 
T„ .in Eq. (11) are . averaged quantities - readers interested in describing 
with ensemble averages can do so for equilibrium as in the literature [53], but 
in this paper the variables are considered to be time averages. 

2.1. Differential form of the diffusive flow 

The diffusive flow J^^ is the difference between the forward and backward 
(one-directional) diffusive flows, e.g., = {^jl^^ - VnO = "^^n^ where the 
arrows indicate migration up or down a potential energy gradient. Using 
Eq. (11), the diffusive flow occurring down a potential energy gradient is 



1 



^T„(r) 
1 



xXn{s) 



Pnir) 

Pn{s) 



m 



m 



(15) 



Equation (15) can be written in differential form by using a truncated Taylor 
expansion to rewrite the variables at site s as variables at site r, 



a„(s) « CLn{r) + A 



dx 



(16) 



where the general character a„ can be any variiable; x, y, and z are the di- 
mensions; m includes n\ whereas the distance X separates the centers of sites r 
and s\ 
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X < 




dx 



kT{r) 



dT 



dx 



y,2 



T„(r) dx 



(17) 



In obtaining this result, nonlinear terms in X were neglected, e.g., the product 
X'^[{doii/dx)y^^{da2/dx)y^^]; this implies that X{dan/dx)y^^ is small but non- 
negligible when compared to a„(r). It should be repeated that the distance A 
unconditionally is constant at a value of unity multiphed by the spatial units 
used in evaluating the derivative; the mean free path should not be equated 
with X because the definition of the spatial derivative in Eq. (16) would be in- 
valid for a varying total density. In this work, we have chosen X to be the dis- 
tance between the centers of sites r and 5. The change of nomenclature also 
should be noted: J^^ becomes J^{r) when writing Eq. (15) in differential form 
- the net flow J^{r) occurs from site r in the x direction, and should not be 
confused with the total flow into a site at position x. Applying Eq. (16) to 
and removing the nonlinear terms in X also gives Eq. (17) - the logic of 
the piecewise Eq. (12) exists in the gradients, which can be positive or nega- 
tive. Hence, the remainder of this paper only discusses the diffusive flow J^{r) 
and the diffusive flux J^{r), i.e., the diffusive flow per area. 



2.2. Existence of Tn and the path dependence of diffusive flux 

Since the relationship between the diffusive flux and the gradient in Eq. (9) is 
linear, the objective of this work is to understand the property r„ and its gra- 
dient. This section constructs the property r„ from the LDFT approach using 
Eqs. (10)-(14), which have been successful at describing diffusion and phase 
behavior as outlined after Eq. (9). Subsequent sections characterize r„ and 
show that it produces the well-known diffusion Eqs. (l)-(4), a vanishing dif- 
fusive flux at the critical point, and agreement with the NET-DPE approach. 

The diffusive flux of species n in the x direction, J^, is one of the components 
ofEq. (9), 



7^ = _ A ^ 
" ~ A dx 



(18) 
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where A = X^ for a cubic lattice. Combining Eqs. (17) and (18) through 



5i; 

where 



dTndT 
dT dx 



d(p^ dx 



dxn dx 



dFn dp„ 



dp„ dx 



y,2 



(19) 



m 



^ = A. / 1 _ V 



1 



^nmPn 



kT ' ^ kT 



m 



m 



1 

— + 



din X.xl 



(20) 



In Eq. (19) and the subsequent equations, partial derivatives indicate the var- 
iables being held constant unless there is no indication, iii which case all var- 
iables are held constant except the two variables in the differentiation. In Eq. 
(20), the subscript '7" represents any species in the set m except for species 
n. 



Equation (19) is a component of the differential of T;,, 



aVn = — 
ox 



x,z 



dz 



(21) 



which is written as dYn because it has a nonexact differential form, i.e., sec- 
ond partial derivatives of r„ are not interchangeable: 



S(Pn d-Cn d(p„ dp„ dT dT dp^ 



dx„ dT ^ dT dT„ 



dpi dp„ dp„ dp, 



(22) 
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This nonexact differential form of Eq. (21) is an important mathematical con- 
sequence of a physical mechanism: diffusion between two points depends on 
the path between these points. Thus, f Jr„ is a path integral, which can be 
evaluated along a particular path of 'avoided^;^j^smg an integrating factor. 

2.3. Integrating factors in non-equilibrium processes 

The classical thermodynamic example of using integrating factors is in Clau- 
sius's theorem, which gives the relation dS = dQrev/T for reversible adiabatic 
processes [54], In this example, dS is the differential change in the system's 
entropy, dQreu is the differential amount of heat added reversibly to the sys- 
tem, and ris the system's temperature that is approximately equal to the heat 
source's temperature. The integrating factor, l/T, modifies dQrev^ creating 
the exact differential dS that can be integrated using only the limits of inte- 
gration. This calculus shows that the heat-transfer process depends on the dif- 
ference in a state variable as well as on the path between the limits that define 
this difference; in particular, the process depends on the path's values of T 
(for an irreversible process, it would be the temperature at which heat is sup- 
plied). Although integration of the resultant exact differential can be advanta- 
geous, an equivalent benefit is the identification of a continuous function 
with continuous partial derivatives that is very relevant to the process under 
investigation. 

In this respect, diffusion is analogous to the heat-transfer process and, there- 
fore, the calculus has similar features. The nonexact form of dFn implies that 
the diffusive flux depends on the path between two limits - integration of dTn 
is impossible unless additional information or integrating factors exist. This 
section shows (i) how diffusion is path-dependent, and (ii) that the driving 
force for the diffusive flux is related to the chemical potential gradient 
through an integrating factor. 

Within the framework of LDFT, an integrating factor exists for dTn^ yielding 
a well-behaved function. Thy that has an exact differential. 



dv, 



n 



(23) 




where the index m includes n. This differential JT„ is 




(24) 



where the fugacity [52] is 
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m 

m "" "■ '"•'-•.vv.. 



and where h = for monomers on a cubic lattice v^tlT lattice spacing h 
Equation (24) is obtained by combining Eqs. (21) and (23) then by integrating 
and using the definition of the fugacity - the partial derivatives in Eq. (21) are 
Eq. (19) and the and z-analogs of Eq. (19). 

Combining Eqs. (9) and (23) gives 

'^'' = ~^^^" = -^^''(l-E/'-)v^«. (26) 

where the concentration C„ = pJX^. Since /„ = P„ exp^^^^^j, it follows that 
VT„ = i^V//„ for an isothermal and isobaric system in which t„ = T„(r); 
thus, the chemical potential gradient appears to drive diffusion. However, 
VT„ also can appear as other gradients. Appendix A shows how Eq. (26) pro- 
duces the well-known diffusion Eqs. (l)-(4). The next section shows that VYjt 
are the conjugate force vectors Xk in the NET-DPE approach. 

2.4. Comparison with NET theory 

The NET Eq. (7) adapted for conservative forces can be rewritten , without 
additional assumptions as 



Q vr 



where fk = Proceeding in the same manner as de Groot and 

Mazur [23] in obtaining Eq. (8), one can amend the definition of the entropy 
flux to 

k 

in order to convert Eq. (27) to 



O" = — 



+ •••, (28) 



where the conjugate force vectors appear in brackets. Near thermodynamic 
equihbrium and in the absence of chemical reactions, it follows that the con- 
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jugate force vectors in Eq. (28). also can be written as 




(29) 



because V In t„ = 0 near equilibrium, V In fcjt = 0 for a constant molecular 
size, and V In 4; = 0 since the coordination nimiber is a constant. A compari- 
son of Eqs. (24) and (29) shows that VTa: = Xr. 

It should be noted that the current LDFT approach, which is for systems that 
are not necessarily isothermal, does not provide equationis for heat flux. Con- 
sequently, it is not yet possible to check whether Onsager's reciprocity rela- 
tions will be obeyed, or whether LDFT correctly predicts the Dufour effect. 
Nevertheless, the analysis up to Eq. (26) also holds for the isothermal case, 
which has the same integrating factor for dTn shown in Eq. (23). The exam- 
ples in the following subsections use the isothermal criterion (as noted) in 
order to provide physical insight to the property r„. 

2,5. Analytical examples of the driving force r„ 

It is useful to know that the function exists with a differential form of Eq. 
(24), with a relationship to the diffusive flux through Eq. (26), and in agree- 
ment with the NET approach with regard to the conjugate force vectors for 
diffusive flux. However, the objectives are to understand the property and 
to provide Eq. (9). This section provides an analysis of F;,. 

Combining Eqs. (23) and (24) gives 



where the latter incorporates the definition of the fugacity, Eq. (25). The fol- 
lowing substitution of variables simplifies this problem and, later, gives more 
insight into the property F;,: 



where 7„ is a dimensiordess ratio representing repulsive to attractive intermo- 
lecular. interactions for species n, and represents the impingement rate of 




(30) 



or 




(31) 
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molecules n onto vacancies at each of the adjacent lattice sites 
substitution, Eqs. (30) and (31) can be rewritten as 

where the ratio X„/ equals fnb/xx„kT. 

We chose this substitution because, for processes of constant 7„, it immedi- 
ately is evident that dTn = dX„ and /„ = -i^X„ where X/A is always con- 
stant; i.e., the diffusive flux is driven by the gradient of the impingement rate 
onto vacancies, at least for constant Y„ processes - this is the case for pure or 
multicomponent systems of interacting molecules under isothermal condi- 
tions with constant V/„ and t„ = r„{T). Even after the substitution, however, 
the second partial derivatives of r„ with respect to X„ and Y„ are not 
interchangeable - integration of dr„ generally requires a "path" between two 
limits. The foUowing sections provide a few reasonable "paths" that charac- 
terize r„. 



Path 1. Near-ideal gases 

For pure near-ideal gases (NIG) and near-ideal gas. mixtures (NIGM), Y„ = 1 
and X„ = -^, because intermolecular mteraction energies are small and the 
insertion probability approaches unity; this also is evident from the high- 
temperature and low-density limit of Eq. (32), which has the insertion proba- 
bility il-J2^pJ. Thus, integration of Eq. (33) gives rj^^rc ^^^igm ^ 
^n+^„> where the term T* is independent of X„. Since there is no diffusion 
of species n when p„ = 0, it is reasonable that the constant r„* = 0, giving 

yNIG _ Y^NIGM _ Pn ^ ( ^2 Pfi \ 

" -"A^vmmr ^^^^ 

which makes use of p„ = P^)? jkT and t„ = Xy/nM„/^kTi where P„ is the 
^MTCAf P^®^^^®' molecular mass of species n. Thus, T^^^ and 

r„ represent the impingement rate of molecu les on their surroundings - 
recall the Hertz-Knudsen formula Pn/y/2nMJcT for the molecular flux to a 
surface in units of \/{area ■ time). Appendix B explains why this section con- 
siders near-ideal gases rather than true ideal gases. 

Path 2. Rigid spheres (one component, isothermal process) 

For a one-component system of rigid spheres under isothermal conditions, 
Y={\-p) and JT = p{\ - p)/xx. Assuming that t = x{T), where M is the 
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molecular mass, allows integration of Eq. (33) to give F^^ = r^^*{T) 
where T^^*{T) is the constant of integration. For all temperatures, however, 
the low-density limit recitiires.i;^f * ^ 0, which is reinforced by the notion that 
the rate of impingement is zero when there are^no^molecules. Consequently, 

XX X 

wh ich make s use of the approximations P{D — b)/kT=l and t = 
a^nM/SkT, where a is the sphere diameter. The term (1 + Pb/kT) in Eq. 
(35) is a correction to the Hertz-Knudsen formulia; for a given pressure, 
fewer molecules impinge on a surface in a period of time for rigid spheres 
than for a near-ideal gas. 



b/a P 

{\+ Pb/kT) y/2nMkf. 



(35) 



Path 3. Rigid sphere mixture (isothermal process at constant total 
density) 

For an isothermal, multicomponent rigid-sphere system at a constant total 
density, Y„ = {I -J2mPm)^ and ^„ =/?„(!- Em /'m)/«T„. Assuming that 
T„ = z„{T) allows integration of Eq. (33) to give F^^ = r^^*{T,p„^„) 
where T^^*{T,p^^„) is the constant of integration that may depend on tem- 
perature and all densities except p„. For all temperatures and densities of 
other species, the low p„ limit requires that r„^^* = 0, which is reinforced 
by the notion that the diffusive flux of species n is zero when there are no mol- 
ecules of species n. Consequently, 

which makes use of th e approximations P{o - b)/kT ^ 1^ == Ylm^rnbmy 
and Xn — Cny/nMn/SkT; the term cOn is the molecular fraction of species n. 
This result represents the impingement rate of species n onto surrounding 
vacancies. 

A major difference among Eqs, (34)-(36) is the appearance of the insertion 
probability (1 - J2mPm) ^ Eq. (36). This is a consequence of the constant 
total-density criterion of Path 3. Equation (36) does not reduce to Eq. (35) in 



bn/(^n 



COnP 



^yinMnkT 



(36) 
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the limit as the number of species decreases to one because Path 3 is restricted 
to a constant total density whereas Path 2 is not. This is explained further in 
Appendix C, and it is an example. Qf how r„ depends on the probabihty of 
there being malleable constituents in thie stirroundi^gs (i;e., vacancies and 
molecules of species whose density gradients can be influenced by diffusion). 



Path 4. Lattice fluid (one component, isothermal process) 

For a one-component lattice fluid with nearest-neighbor interactions and no 
external fields, Y = (1 - pf exp{-xep/kT) and X = p{l -p)lx'c. Assuming 
that T = x{T) allows integration of Eq. (33), 



r = -!- 




(37) 



which makes use of the limit T{p = 0) = 0. When there are no attractions 
(e = 0), Eq. (37) gives the one-component rigid-sphere result p/xx in Eq. 
(35). This result does not have the insertion probability that appears in the 
numerator of Eq. (36) because this path does not have a constant total- 
density criterion. The lack of an insertion probability is explained further in 
Appendix C and it is an example of how r„ is not a simple function. 

As a test of Eq. (37), combination with Eq. (9) under the restrictions of this 
"path" gives 



7 = -l 

XT 



kT 



VC, (38) 



where C — pf)? and the area through which diffusive flow occurs is = A^. 
This is the proposed diffusion equation that gives figures 2-6 in reference [50]. 
When T =Tc and /? = i, the bracketed term in Eq. (38) equals zero - this is 
the correct result at the critical point [2]. Thus, the new driving force provided 
in Eq. (37) gives a vanishing diffusive flux at the critical point, i.e., it gives 
item (/) as promised in the introduction. Furthermore, the bracketed term 
has negative values for states below the spinodal function on a T - p phase 
diagram; i.e., diffusion occurs against concentration gradients in the spinodal 
decomposition process. 

The driving force in Eq. (37) also gives the fluxes and non-equilibrium phase 
behavior in figures 2-6 of reference [50], which yields the non-equilibrium 
density profiles in the x dimension by solving a steady-state material balance 
equation to give the profile 
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X 

L 



kT\2 



kT\2 



Pl 



It now is evident that this profile from reference [50] is an expression for 
the Unearity of the driving force in Eq. (37); i.e., it is equivalent to 

rr(L) - r(o)i 



r{x) = 



x + r{0). 



The discussion in this paper provides a very simple physical interpretation of 
the driving force in Eq. (37). 

2.6. Graphical analysis of the driving force T„ 

Equations (34)-(36) show that r„ represents a rate of impingement of mole- 
cules onto malleable constituents in the surroundings, at least for near-ideal 
gases, rigid spheres, and constant total-density rigid-sphere mixtures. This is 
not obvious for the lattice fluid undergoing Path 4, i.e., Eq. (37) appears to 
be an impingement rate only in the low-density limit. For a different perspec- 
tive of r„, this section provides a graphical analysis by using Eq. (33), which 
shows that r„ is the area under the "path" in Y„ versus X„/Y„ space. Recall 
from Eq. (32) that Y„ is a ratio of intermolecular repulsions to attractions and 
that X„/ Y„ = fnb/zr„kT. 

Figure 1 shows examples of Paths 1-4 in Y„ versus X„/Y„ space. At low 
values of X„/Y„, the area r„^ > > > In going from Path 1 to 2, 




4 & 8 10 12 

Figure 1 Isotherms in versus X„/Yn space. Paths: (1) pure near-ideal gas and near-ideal 
gas mixture, (2) rigid spheres, (3) rigid sphere mixture at constant total density, and (4) lattice 
fluid. The foUowing parameters were chosen arbitrarily: x„ — 8.5 x 10"^^^ and = 6 for all 
paths; YlmPm = for Path 3; and e/kT = 0.5 for Path. 4. Areas imder the curves equal the 
property r„, which is the rate of impingement of molecxiles on their surroundings at least for 
Paths 1-3 as described in Eqs. (33)-(36). 
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the decrease in area is caused by the rigidTsphere repulsion of species n; in 
going from Path 2 to 3, the decrease is caused by the rigid-sphere repulsion 
from all the species in the mixture-- the area under Path 3 is 65% of the area 
under Path 1 due to the insertion probabiUty 0:65.= (1 - J^^p^). On the 
other hand, in going from Path 1 to 4, the increase in area at low Z„/ 7„ is 
caused by a greater contribution from attractive interactions than from repul- 
sive interactions, which dominate at higher values of X„/Yn for the lattice 
fluid. If this low X„/Yn behavior appears counterintuitive, consider that in 
the low X„/ Yn range this lattice fluid has a greater pressure than an ideal gas 
with the same fugacity; i.e., (f/P) < 1 for the lattice fluid at low density, 
whereas / = P for an ideal gas. Therefore, at low X„/ Y„, the rate of impinge- 
ment of the lattice fluid is greater than that of the ideal gas. 

Path 4 is the only path in Figure 1 that does not have monotonic behavior 
and it is the "path" used in reference [50]; therefore, it is worth investigating 
it further. Figure 2 shows more examples of lattice fluids undergoing Path 4. 
For fluids that are below their critical temperatures, the path can be multi- 
valued at a particular value of X„/ Y„ - this is not the case for fluids that are 
above their critical temperatures. Figure 3 describes the meaning of this non- 
trivial behavior of subcritical fluids undergoing Path 4 for a particular value 
of e/kT. 

For interpreting Figure 3, it may be useful to consider an isothermal lattice 
flmd in a piston-cylinder system. At vapor-Uquid equilibrium (VLB) there is 
no diflfusive flux; this occurs at the intersection points 2 and 8 in Figure 3, 




1 2 3 4 5 



6 



Figure 2 Lattice fluids with different values of e/kT undergoing Path 4 in Y„ versus Xn/Y„ 
space. For all data, r„ = 8.5 x lO'^s and « = 6. The bolded curve (f) is the Path 4 curve 
shown in Figure 1. The fluid giving curve (d) has an e/kT = -2/3, which represents the critical 
temperature; the maximum along this curve occurs exacdy at the critical point. The fluids giv- 
ing curves (a)-(c) are below their critical temperatures, whereas the fluids giving curves (e)-(h) 
are above theu: critical temperatures. 
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4,0 



3.0 



2.0 



1.0 





1 2 .3 4 5 

Figure 3 Subcritical lattice fluid .undergoing Path 4 in y„ versus Xn/Yn space. For this path, 
x„ = 8.5 X 10"^^j, e/kT = -0.88, and = 6; it is exactly the same as the corresponding curve 
in Figure 2. Subplots show the locations of points on the phase diagrams. In the kT/\£\ versus 
p diagram, the solid curve is the binodal and the dashed curve is the spinodal; the critical point 
is at kT/\e\ = 1.5 and p ~ 0.5. Points 2 and 8 are at the same location in Yn versus Xn/Y„ 
. space, but not on the phase diagrams. 




1.0 1.2 1.4 1.6 1.8 



X„/Y„ 

Figure 4 Dlustration of the differences r„"^(7) - F^pO) and r7P(3) - F„"^(9) for the lattice 
fluid described in Figure 3. When the system is pressurized to the metastable state represented 
by points 3 and 9, diffusive flux occurs from the vapor to the Uquid phase; for the metastable 
state represented by points 1 and 7, diffusive flux occurs from the hquid to the vapor. When 
integrating from point 1 to 7, the area mostly is positive (-h) except for the small region labeled 
negative (— ); ld the integration from point 8 to point 7, d(Xn/ Yn) < 0. 

which shows that the areas under the "path" to these points are the same, i.e., 
the driving forces r^^^^(2) = r^'^(8). When this system is pressurized sUghtly 
to the metastable state represented by points 3 and 9, diffusive flux occurs 
from the vapor to the Uquid phase because r„''^^(3) > r„^^(9); this reduces 
the pressure, returning the system to VLE. On the other hand, if the system 
originally in VLE is brought to the metastable state represented by points 1 
and 7, a diffusive flux would occur from the Uquid phase to the vapor phase 
because r„^^(7) > r;^P(l). Figure 4 shades both differences r;^P(3) - r^^^(9) 
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1.0 ■ 1..2 • 1.4 . 1.6 1.8 



Figure 5 Illustration of the difference r„'^'*(9) - T^^Cl) > 0; for the lattice fluid described in 
Figure 3. This fluid encoxmters a vapor-liquid phase transition, as described m reference [50], 
in. which diffusive flux occurs from the liquid to the vapor phase. Points 1 and 9 are the bound- 
ary conditions, which are in the stable one-phase region of the phase diagram, 

and r^^^*^(7) - r„''^^(l); note that integration from point 8 toward the points 7 
arid 6 produces negative areas, i.e., r^''*{8) > r!^\7)> r„^'*^(6). 

In figure 6 of reference [50], two isothermal diffusion processes encountering a 
vapor-liquid phase transition are summarized: diffusive flux occurs (i) from 
the Uquid to the vapor phase when the boundary conditions represent fluids 
in the stable one-phase region, and (ii) from the vapor to the hquid phase 
when the boimdary conditions represent fluids in the metastable regions. An 
example of case (i) is a two phase system at points 1 and 9 in Figure 3; diflfu- 
sive flux occurs from the Uquid to the vapor because r^^(9) > r^^^{l) as 
shaded in Figure 5. On the other hand, an example of case (ii) is a two-phase 
system at points 3 and 7 in Figure 3; diffusive flux occurs from the vapor to 
the hquid in this case because r„^^^(3) > r„^^(7) as shaded in Figure 6. 

As a final example, consider the phenomenon of spinodal decomposition. 
When a uniform supercritical fluid is quenched to temperature below the spi- 
nodal function in the subplot of Figure 3, the system immediately becomes 
unstable; a small fluctuation in p due to thermal motion causes decomposi- 
tion of this system into two stable phases. Figure 7 illustrates the driving 
force behind this phenomenon for a lattice flxiid originally uniform at /? = 0.5 
and quenched to kT/\e\ = 0.714 (the critical point occurs at kT/\e\ = 1.5). In 
this figure, Ap is a small fluctuation in the density, which can be caused by 
thermal motion and which creates the driving force difference r{p — Ap) — 
T{p + Ap). Note that integration along the path in the negative Xn/Y^ direc- 
tion gives a negative value for the area; therefore, diffusion occurs from 
regions of lower density to regions of higher density because r{p — Ap) > 
r{p + Ap). This enlarges the density fluctuation Ap and, consequently, the 
difference r{p - Ap) - r{p + Ap), accelerating the decomposition. 
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Figure 6 Illustration of the differences r^^^{3) - T^'^il) > 0; for the lattice fluid described in 
Figure 3. This fluid encounters a vapor-hquid phase transition, as described in reference [50], 
in which diffusive flux occurs from the vapor to the Uquid phase. Points 3 and 7 are the bound- 
ary conditions, which are in the nietastable regions of the phase diagram. 
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Figure 7 Driving force for spinodal decomposition of a lattice fluid originally uniform at 
p = 0.5 that is quenched from supercritical temperature to a temperature kT/\sl= 0.714 that 
is below the spinodal function (the critical temperature is kT/\e\ = 1.5). The driving force is 
the area r„{p - Ap) - r„(/? + Ap) > 0, where A/? is a small fluctuation in the density; diffu- 
sion occurs from lower to higher densities, enlarging the fluctuation and accelerating the 
decomposition. 

3. Discussion 

Numerous equations have been used to model diffusive flux, but none of these 
equations generally have a linear relationship between the diflftisive flux and 
the gradient; the coeflicients generally depend on the state variables, which 
usually depend on position and time. This is a mathematical consequence of 
, diffusion's path-dependence. The diffusive flux depends on more than just the 
difference in the state variables, much like how the heat flux in Clausius's the- 
orem depends on the integrating factor 1/r in addition to the difference in 
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the entropy. Since diffusion is analogous to heat transfer, an integrating fac- 
tor exists and relates the diffusive flux to a difference in state variables (e.g., 
to a difference in chemical potentiaJU)r fugacity). In the LDFT framework, 
this is summarized by Eq. (26), wMcK'uses the.in^grating factor that multi- 
pUes dr„ in Eq. (23); the chemical potential gradienT^Ould arise from the 
gradient in Eq. (26) by imposing the restrictions of an isothermal and isobaric 
process in which t„ = t„(7'). 

If one integrating factor exists, however, there exist an infinite number of 
integrating factors that can have various forms. For example, consider 
Eq. (31), which shows that e(?'"+«'')/*^^/(l - EmP/i.)^ also is a valid integrat- 
ing factor in addition to that in Eq. (30). This highlights a fundamental 
problem for diffusion: the driving force for the diffusive flux can appear to 
be a variety of gradients. For example, both V/„ and V In /, can appear to be 
driving forces; the latter is equivalent to ^.VMh for an isothermal and iso- 
baric system. Although the conjugate gradient may change due to the choice 
of integrating factor as demonstrated by Eqs. (!30) and (31), the gradient 
Vr„ remains independent of the choice of integrating factors. Furthermore, 
Vr„ always is linearly proportional to the diffusive flux, as shown in Eq. 
(9). Therefore, we propose that Vr„ is a more general thermodynamic driv- 
ing force for diffusion. The challenge is determining what r„ represents 
physically. 

Equations (34)-(36) and Appendix C show that r„ appears as the impinge- 
ment rate of n molecules onto the malleable constituents in the surroimd- 
ings; at least for near-ideal gases, rigid spheres, and constant total-density 
rigid-sphere mixtures. This is reassuring because kinetic theory approaches 
for gas diffusion begin with defining the flux as the difference in the rate of 
impingement [55]; these approaches are for constant total-density systems 
in the dilute limit, which has an insertion probability of unity - the proba- 
bility of encountering a malleable constituent is unity. For the lattice fluid 
undergoing Path 4, however, the impingement rate is not obvious in Eq. 

(37) outside of the low-density limit. This is not necessarily a mistake - 
the partial derivative of T with respect to p gives the diffusive flux in Eq. 

(38) , predictmg the correct (vanishing flux) behavior at the critical point and 
the non-equilibrium phase transitions in reference [50]. Of course, this mean- 
field result is qualitative, and it would need to be extended for quantitative 
accuracy [47]. 

The property r„ can be difficuh to interpret because it is not a state function 
and it is not apparent unless a "path" is specified. The following is one very 
simple physical interpretation of t that clarifies why the impiiigement rate is 
not obvious in Eq. (37) outside of the low-density limit. First consider the 
one-dimensional case of Eq. (38), 
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^ Ik" (kT .jjmi.y_c 

^ ~ ^xkT\l ) dx J 

where At/'"''^ = -xep{\ - p)/2.and is'tiie' ehan^j^ energy 
per site upon mixing for regular solutions of vacancies and molecules with 
nearest-neighbor interactions, and the term kT/l is the x component of the 
mean kinetic energy per molecule. The term HU^P^ is positive for intermolec- 
ular attractions, acting as an energetic penalty for the mixing of molecules and 
vacancies, and contributing to vapor-Uquid phase transitions at low enough 
temperatures. Therefore, when molecules' kinetic energy in the direction ex- 
ceeds the energetic penalty for interchanging with vacancies, diffusion occurs 
down the concentration gradient in the x direction (i.e., molecules move to- 
ward the vacancies). If AC/^"'^ is greater in magnitude than the mean kinetic 
energy, diffusion occurs against this concentration gradient. To make the 
connection to the property T, consider the case kTIl > AU^^. When non- 
interacting molecules impinge on an imaginary surface perpendicular to the di- 
rection of motion, on average they carry kT/2 through the surface in a short 
period of time. When intermolecular attraction occurs, the amount of inter- 
nal energy carried through this imaginary surface is less than kT/2, as there 
is a toll AU^^ for moving into a vacant site. The dimensionless amount of 
intemal energy per site that passes through the surface in a period of time is 

Thus, r can be interpreted as a rate of dimensionless intemal energy that is 
available per site for molecular translation in each of x directions. The prod- 
uct rkT appears to be the available power for molecular translation in the 
±x directions. For near-ideal gas molecules and noninteracting rigid spheres 
{AU^^ = 0), this rate of dimensionless energy appears as the rate of impinge- 
ment in each of x directions: 

^^p kT/2 ^ p 

XT kT/2 xr' ' • , - 

This is one interpretation of Path 4. 

There can be an infinite number of paths between two thermodynamic states. 
Many isothermal paths involving steady-state diffusion of interacting species 
have constant fugacity gradients; even processes involving mixtures of inter- 
acting molecules in a strong extemal field exhibit nearly hnear fugacity pro- 
files in the dimension of the diffusive flux [52]. Such processes have constant 
values of in Eqs. (32) and (33). For. these cases dFn = dXn, i.e., the diffu- 
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Figure 8 Comparison of F, fx, and / for the lattice fluid with s/kT = -0.88 and = 6; all are 
shown in dimensionless form versus the dimensionless pressure. A dimensionless pressxire of 
0.100 is equivalent to 10.5 MPa for argon at a temperature of 300 K. Values of the dimension- 
less chemical potential appear on the right axis in excess of the ideal gas reference potential 
Values of the dimensionless T and / appear on the left axis. 

sive flux is driven by the gradient of the impingement rate onto vacancies. 
Figure 8 compares r„, ju^, and for a lattice fluid undergoing Path 4. These 
properties exhibit very similar behavior; therefore, it is not surprising that r„ 
can be mistaken for an energy density such as the chemical potential or the 
fugacity. We hope that the skeptical reader would at least recognize that Eqs. 
(3) and (4) are identical and that the coefficient to V/„ is constant more often 
than the coefficient to V//„. This occurs mostly because it is more probable to 
have a uniform Qot than a uniform C„ when there is diffusive flux of species «. 

Equation (9) is a linear law describing diffusion; the linear relationship be- 
tween the diffusive flux and tke gradient remains in the presence of gradients 
'^^Pny Pm^^ <Pn^ Pm fn, T, and T„. The gradient in Eq. (9) transforms into a well- 
recognized gradient (e.g., the chemical potential gradient) only when par- 
ticular properties are uniform, as shown in Appendix A. If a sufficient num- 
ber of properties are uniform, the resultant diffusion equation can have a 
linear diffusive flux with respect to any one of several different gradients. 

4. Conclusions 

The LDFT framework was used to formulate a linear diffusion equation that 
has a constant coefficient even in the presence of gradients in all the state 
variables. It is proposed that the gradient in this equation is the driving force 
for diffusion. This driving force is related to the gradient of state variables 
(e.g., to the chemical potential gradient) through integrating factors, much 
like in Clausius's theorem, which relates the heat flux to the difference in en- 
tropy through the integrating factor l/T. At constant temperature, for near- 
ideal gases, rigid spheres, and constant tofal-density rigid-sphere mixtures, 
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the driving force for diffusion appears to be the gradient of the impingement 
rate of molecules on malleable constituents of the surroundings. For the lat- 
tice fluid with nearest-neighbor attractions, the driving force appears to be the 
gradient of the impingement rate at'ldw densiti^; but it deviates outside of 
this limit, appearing more like a rate of dimensioniess internal energy avail- 
able for molecular translation - this rate of dimensioniess internal energy 
equals the rate of impingement when the configurational energy of mixing 
molecules and vacancies is zero. This paper shows that the proposed driving 
force gives the well-known diffusion equations, a vanishing diffusive flux at 
the critical point, and agreement with the theory of non-equilibrium thermo- 
dynamics. Furthermore, this paper uses this driving force to explore various 
diffusion equations, finding that the diffusive flux is linear more often with the 
fugacity gradient than with the chemical potential gradient because, when 
species n diffuses, it is more probable to have a constant total-density than a 
constant density of species n. 

Appendix A: Unifonn properties produce the well-known diffusion 
equations 

Equations (l)-(4) are simplifications of a more general equation for the diffu- 
sive flux. Fickian behavior (with a constant diffusivity) is observed only when 
a sufi&cient number of state variables are uniform. This section demonstrates 
this by showing that Eq. (26) reduces to the well-known diffusion equations. 

For isothermal and isobaric systems m which t„ = T„(r), Eq. (26) becomes 

which is Eq. (4) in the LDFT framework. On the other hand, for isothermal 
systems with t„ = T„(r), Eq. (26) has the form of 

which is Eq. (3). Equation (A.2) simpUfies further when the fugacity coeffi- 
cient ^„ is constant, 

which is Eq. (2). Lastly, for isothermal systems with t„ = t„(T) and the addi- 
tional restrictions of a constant total density and of no gradients in (p„ and £«, 
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Eq. (26) has the form of 

7, = -iw„ = - --&.^y(^,_^ (A.4) 

which is Eq. (1). 

Again, it should be noted that these results are for lattice fluids. To make 
these results more general, the term (1 - ^„pj should be replaced with the 
insertion probability as shown in reference [50]. 



Appendix B: Near-ideal gases versus ideal gases 

Path 1 describes near-ideal gases, i.e., real gases that experience collisions but 
which have properties approaching ideal gas behavior characterized by the 
equation PV = nRT. On the other hand, ideal gases are point masses that 
do not experience collisions and, therefore, pure ideal gases cannot have an 
average velocity that differs from the mass-averaged velocity; i.e., their diffu- 
sion velocity [56] is zero at uniform temperature and mass-averaged velocity: 

F, = (D,-u,)=ljFAfN^-, (B.l) 

where the diffiision velocity F,- is the average velocity of i-molecules y,- in 
excess of the mass-averaged velocity Vo, Vt is the particular velocity of an 
i-molecule with respect to the mass-averaged velocity, m is the number of 
z-molecules per' volume, the (Maxwell-Boltzmann) distribution ff^ is the first 
term in the velocity distribution function in the Boltzmann equation [57], and 
<t>i is the perturbation function that gives the first-order perturbation from the 
Maxwell-Boltzmann distribution. For gases that experience coUisions, can 
be evaluated and is related to the coUision integral that appears m the diffu- 
sivity for pure or binary gases at low density [55]. On the other hand, for an 
isothermal system of noninteracting point masses with a uniform mass- 
averaged velocity, 0)/^ = 0 giving Kf^ = 0. For this reason, Path 1 considers 
near-ideal gases, which undergo collisions and have a velocity distribution 
and, therefore, can have diffusion velocity. 

Appendix C: "]VIalleabIe constituents" and the insertion 
probability in r„ 

The property r„ has different forms that depend on the constraints/conditions 
of the diffusion "path". For example, the insertion probability (1 - J^^Pm) 
seems to appear arbitrarily in Eq. (36) but it does not appear in Eqs.lsi), 

J. Non-Equilib. Therinodyn. • 2006 • Vol. 3 1 • No. 4 



380 



D. Matuszak et al. 



(35), and (37). This section demonstrates that the foms of have a common 
origin; that the reason why the insertion probabihty appears in Eq. (36) is the 
same for the lack of an insertion probabihty in Eqs, (34), (35), and (37); and 
that is not limited to the foms Of Eqs. .(3|)^^^ 

Previously, it was stated that r„ can appear as the impingement rate onto the 
malleable constituents in the surroundings. For lack of a better term, we 
consider malleable constituents to include vacancies and molecules of species 
whose density gradients can be influenced by diffusion; e.g., a constant- 
density matrix of a solid membrane does not qualify as a malleable constitu- 
ent even if the molecules of this matrix wander; this is because the density of 
the matrix does not respond to penetrants. Similarly, a multicomponent sys- 
tem that is restricted to a constant total density, overall, is considered a ma- 
trix that does not respond to individual penetrants within the system. With 
this as the context, we provide the following illustrative example. 

Consider a binary system of identical nonattracting spheres that differ only 
in their "color". Spheres are labeled ''1" and "2" as in the problem of iso- 
thermal color counterdiflfusion of Weeks-Chandler-Anderson fluids [33, 52]. 
The density profiles of both species are Unear in the absence of external fields, 
/?„ = mnX + bn where wi^ is the slope and bn is the mtercept. For this problem 
we have 



There are at least three interesting cases of Eq. (C.l). First, Fj^ — 
^{i - Pi - P2) for equal color coimterdiffusion at a constant total density 
(m2 == -mi); this is the binary version of the constant total-density rigid- 
sphere mixture result of Eq. (36) - the probability of having a malleable con- 
stituent is the vacancy probabihty p^=zi -p^- Second, F/^ = ^ (1 - p^) 
for diffusion of species 1 through a matrix of species 2 in which V/?2 = ^2 = 
0; note that the insertion probability does not appear explicitly - the proba- 
bihty of having a malleable constituent is the sum = 1 — /?2- Third, 
r/^^ = ^ for diffusion of a pure fluid in which half the molecules arbitrarily 
are labeled "1" {mi = m2 and p^ = P2); i.e., F^^^ = ^ where p = 2p^ for the 
pure fluid; this is the rigid sphere result, Eq. (35), as well as the noninteracting 
lattice fluid result, as can be seen from Eq. (37) - the probability of having a 
malleable constituent is the sum p^+p = I, 



= x{T), Y„ = (1 -p, -P2)\ X„ = ^{l-pi -p,). 



This allows integration of Eq. (33) for an isothermal path, 




(CI) 
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Thus, the insertion probability may appear expUcitly in r„ if it equals the 
probability of having a malleable constituent. However, if the insertion prob- 
ability does not appear explicitly,,4t^contributes to the flux and driving flux 
through the equation of state. To reijffof ce' this 4asjtj)o consider that Eqs. 
(34) and (35) give 

r^^^(r,p)>r^^(r,p). 
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Use of L-arginine, L-omithine, or L-citruUine and topical preparations with these substances 

The present invention relates to the use of L-arginine, L-omithine, and/or L-citrulline and the use 
of cosmetic and dermatological topical preparations that contain these active ingredients. 

The epidermis is rich in nerves and nerve end apparatus such as Vater-Pacini corpuscles, Merkel- 
neurite complexes, and free nerve endings for feeling pain, cold, heat, and itching. 

In people who have sensitive skin or skin that is prone to injury it is therefore possible to observe 
a neurosensory phenomenon called "stinging" (in English, "sting" means to injure, bum, cause 
pain). This "sensitive skin" is fundamentally different from "dry skin" that has a thickened and 
hardened stratum comeum. 
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Typical reactions when there is "stinging" on sensitive skin are reddening, tension, and burning 
of the skin, as well as itching. 

Itching in the case of atopic skin and itching in the case of skin diseases is also considered a 
neurosensory phenomenon. 

The object of the invention is therefore to provide active substances, and topical preparations 
having such active substances, that prevent neurosensory phenomena or alleviate such 
phenomena or provide rapid relief for them after they occur, that is, that are suitable for 
prophylactic purposes and/or treatment. 

"Stinging" phenomena can be considered disorders to be treated cosmetically. In contrast, 
severe itching, in particular severe itching of the skin that occurs with atopy, can also be called a 
serious dermatological disturbance. 

These objects are attained by the use of one or a plurality of compounds selected from the group 
of L-arginine, L-ornithine, and L-citruUine or their salts, acid addition salts, esters, or amides, 
where necessary with the addition of folic acid or its salts and/or one or a plurality of compounds 
selected from the group of flavins, for prophylactic purposes and/or treatment of neurosensory 
phenomena. 

They are preferably used for prophylactic purposes and/or treating neurosensory phenomena of 
the skin, in particular "stinging" and atopy (neurodermitis). 

The subject-matter of the invention is also the use of cosmetic and dermatological topical 
preparations that contain the inventive active substances, 
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for prophylactic purposes and/or for treating the neurosensory phenomena. 

One skilled in the art would not be able to foresee that the inventive active substances and that 
cosmetic or dermatological preparations having active ingredients of arginine, citrulline, 
ornithine, folic acid, and flavins in accordance with the invention would be excellently suited for 
prophylactic purposes and treatment of neurosensory phenomena. 

Preferred salts of L-arginine, L-omithine, and L-citrulline are water soluble salts, e.g. sodium, 
potassium, and ammonia salts. This also applies to the acid addition salts. Suitable acid 
additional salts are obtained with inorganic and organic acids. Hydrochlorides, sulfates, acetates, 
caprylates, and citrates are preferred. 

Suitable esters of these compounds are e.g. those that are formed with short-chain and medium- 
chain alcohols, preferably mono-alcohols, but in particular methanol, ethanol, or propanol. The 
ethylesters are preferred. 

Preferred amides are short-chain and medium-chain mono-alkylamides and di-alkylamides. 

Alkyls with the aforesaid substituents contain e.g. up to 12 carbon atoms, preferably up to 6 
carbon atoms. 

Particularly preferred are active substance combinations and topical preparations that contain L- 
arginine and/or its inventive derivatives. 

L-arginine and its derivatives are also distinguished by particularly good skin penetration ability. 
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The inventive amino acids and/or their derivatives are preferably contained in the inventive 
cosmetic and dermatological preparations in quantities of 0.01 to 30 vvt.%, particularly preferred 
0.01 to 10 wt.%, especially 0.1 - 7.5 wt.%, relative to the total preparation. L-arginine, 
ornithine, and citrulline and their derivatives can be employed individually or in combination; 
however in this case the preparations then contain particularly preferred L-arginine in quantities 
of 1 to 10 wt.%, relative to the total preparation. 

Particularly suitable salts of folic acids are water-soluble salts, especially sodium, potassium, and 
ammonia salts. 

Preferred compounds from the group of flavins are flavin adenine dinucleotide (FAD) and flavin 
mononucleotide (FMN). 

Folic acid or its salts and/or flavins are advantageously contained in the inventive preparations, 
preferably each in quantities of 0.0001 to 5 wt.%, especially 0.01 to 1.5 wt.%, in each case 
relative to the total weight of the preparation. 

The inventive cosmetic or dermatological topical preparations can be based on conventional 
formulation basics per se and can treat the skin in the sense of a dermatological treatment or a 
treatment in the sense of skin-care cosmetics. 

Particularly advantageous preparations are furthermore obtained when the inventive active 
substances are combined with antioxidants. 
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It was not foreseeable that the active substances and combinations with an active content of 
arginine, ornithine, or citrulline and the aforesaid co-factors or the aforesaid combination with 
antioxidants would alleviate neurosensory hypersensitivity and itching. Furthermore, it was not 
foreseeable that they would lead to products that are tolerated by skin or that increase tolerance 
and in healthy skin does not interfere with autologous microflora in the skin. 

Thus the objects cited in the foregoing are attained. 

Typical troublesome neurosensory phenomena that are associated with the terms "stinging" and 
"sensitive skin" are skin reddening, pins and needles, pricking, tension, burning of the skin, and 
itching. They can be evoked by stimulating environmental conditions, e.g. massage, the effect of 
surfactants, weather factors such as the sun, cold, and dryness, but also damp heat, heat radiation, 
and UV radiation, e.g. from the sun. 

Surprisingly, in accordance with the invention the symptoms of sensitive skin and itching are 
alleviated or prevented in diseases of the skin such as e.g. atopy. 

The inventive antioxidants can advantageously be selected from the group of usual cosmetic and 
dermatological antioxidants, especially from the group comprising tocopherols and their 
derivatives, especially a-tocopherol or a-tocopheryl esters, in particular a-tocopheryl acetate, 
furthermore sesamol, bile acid derivatives such as methyl, ethyl, propyl, amyl, butyl, and lauric 
gallate, the coniferylbenzoate of benzoin resin, nordihydroguaiacinic acid, nordihydroguaiaretic 
acid. 
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butylhydroxyanisol, butylhydroxytoluene, ascorbic acid, citric acid, phosphoric acid, lecithin, 
trihydroxybutyrophenon, carotenes, vitamin A and its derivatives, especially retinylpalmitate, 
ascorbic acid, ascorbylpalmitate, dilaurylthiodipropionate, distearylthiodipropionate, 
monoisopropylcitrate, thiodipropionic acid, EDTA and EDTA derivatives, cysteine, glutathione 
and esters, uric acid, lipoic acid and esters, carotenes, heavy metal complex formers such as 
delta-amino lavulinic acid, and phytin acid. 

The inventive cosmetic or dermatological preparations contain preferably 0.01 to 10 \vt,%, but 
especially 0.1 to 6 wt.%, relative to the total weight of the preparations, of one or a plurality of 
the substances from the group of antioxidants. 

It is preferred to select the inventive antioxidants from the group of tocopherols and their 
derivatives. 

Moreover, the preparations advantageously also contain urea. The urea content is for instance 
0.01 - 30 wt.%, especially 0.1-10 wt.%, relative to the total weight of the preparations. 

For application, a sufficient quantity of the preparation is applied to the skin in the manner 
normal for cosmetic and dermatological products. 

Especially preferred are skin care and hand care preparations, sun screen and after-sun 
preparations, and bath and shower preparations having a skin-care fiinction. 
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Dermatological and cosmetic preparations in accordance with the invention can come in different 
forms. Thus e.g. aqueous, alcoholic, or aqueous-alcoholic solutions, oil-in-water (0/W) 
emulsions, water-in-oil (W/0) emulsions, multiple emulsions (e.g. water-in-oil-in-water 
(W/O/W) emulsions), gels, hydrodispersions, solid blocks and aerosols can contain the aforesaid 
combinations of active substances. 

The inventive topical preparations can contain the usual adjuvants such as emulsifiers and 
preservatives. 

Those cosmetic and dermatological preparations that are in the form of sun screens are also 
preferred. These advantageously also contain at least one UVA filter and/or at least one UVB 
filter and/or at least one in organic pigment. 

However, those preparations that are applied after the skin has been exposed to light, that is, 
after-sun products, are also extremely advantageous. With such preparations, whether additional 
UV filter substances should be used or not is left to the discretion of the one skilled in the art. 

Cosmetic preparations in accordance with the invention for protecting the skin from UV rays can 
be in different forms, such as they are e.g. normally used for this type of preparation. Thus they 
can be e.g. an aqueous, alcoholic, or aqueous alcoholic solution, a water-in-oil (W/0) emulsion 
or an oil-in-water (0/W) emulsion or a multiple emulsions, for instance water-in-oil-in-water 
(W/O/W) emulsion, 
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a gel, a hydrodispersion, a solid block, or an aerosol. 

The inventive topical preparations can contain cosmetic adjuvants such as are normally used in 
such preparations, e.g. preservatives, bactericides, fragrances, anti-foaming agents, dyes, 
pigments that have a dying action, thickening agents, surfactant substances, emulsifiers, 
softening substances, wetting and/or moisture-retaining substances, fats, oils, waxes, or other 
conventional constituents of a cosmetic formulation such as alcohols, polyols, polymers, foam 
stabilizers, electrolytes, organic solvents, and silicone derivatives. 

If the cosmetic or dermatological preparation is a solution or lotion, the following can be used 
for solvents: 

water or aqueous solutions; 

oils, such as triglycerides of capric acid or caprylic acid, but preferably also castor oil; 

Fats, waxes, and other natural and synthetic fatty bodies, preferably esters of fatty acids 
with alcohols having a low C number, e.g. with isopropanol, propylene glycol, or 
glycerol, or esters of fatty alcohols with alkane acids having a low C number or with fatty 
acids; 

Alcohols, diols, or polyols having a low C number, and their ethers, preferably ethanol, 
isopropanol, propylene glycol, glycerol, ethylene glycol, ethylene glycolmonoethyl or 
monobutyl ether, propylene glycol monomethyl. 
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Monoethyl or monobutyl ether, diethylene glycol monomethyl or monoethyl ether 
and similar products. 

In particular mixtures of the aforesaid solvents are used. Water can be another constituent in 
alcoholic solvents. 

Emulsions in accordance with the invention, e.g. in the form of a sun-screen cream, sun-screen 
lotion, or suntan lotion are advantageous and contain e.g. the aforesaid fats, oils, waxes, and 
other fatty bodies, as well as water and an emulsifier as it is normally used for such a 
formulation. 

Cosmetic and dermatological preparations for treating and caring for the skin can be gels that, in 
addition to the active substances and solvents normally used for them, also contain organic 
thickening agents, e.g. gum arable, xanthan gum, sodium alginate, cellulose derivatives, 
preferably methyl cellulose, hydroxymethyl cellulose, hydroxyethyl cellulose, hydroxypropyl 
cellulose, hydroxypropylmethyl cellulose, and inorganic thickening agents, e.g. aluminum 
silicates such as for instance bentonites, or a mixture of polytheylene glycol and polyethylene 
glycol stearate or distearate. The thickening agent is contained in the gel e.g. in a quantity 
between 0.1 and 30 wt.%, preferably between 0.5 and 15 wt.%. 

Inventive gels normally contain alcohols having a low C number, e.g. ethanol, isopropanol, 1.2- 
propandiol, glycerol, and water or one of the aforesaid oils in the presence of a thickening agent, 
that 
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in oily-alcoholic gels is preferably silicon dioxide or an aluminum silicate, and in aqueous- 
alcoholic or alcoholic gels is preferably a polyacrylate. 

Hydrodispersons are dispersions of a liquid, semi-solid, or solid internal (discontinuous) lipid 
phase in an external aqueous (continuous) phase. 

In contrast to 0/W emulsions, which are distinguished by a similar phase arrangement, 
hydrodispersions are largely free of emulsifiers. Like emulsions, hydrodispersions are 
metastable systems and tend to convert to a condition of two discrete cohesive phases. In 
emulsions, the selection of a suitable emulsifier inhibits phase separation. 

In hydrodispersions of a liquid lipid phase in an exterior aqueous phase, the stability of such a 
system can be assured for instance in that a gel structure in which the lipid droplets are 
suspended stably is constructed in the aqueous phase. 

Solid bars in accordance with the invention can contain e.g. natural or synthetic waxes, fatty 
alcohols, or fatty acid esters. Lipstick-style bars are preferred. 

Suitable as propellants for inventive cosmetic or dermatological preparations that can be sprayed 
from aerosol containers are the normal, known, highly volatile, liquefied propellants, for instance 
hydrocarbons (propane, butane, isobutane), which can be used 
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alone or in mixtures. Using compressed air is also advantageous. 

Naturally one skilled in the art also knows that per se there are non-toxic propellant gases that 
fundamentally would be suitable for the present invention, but that should not be used because 
they may be hazardous to the environment or because of some other factor, especially 
fluorocarbons and hydrochlorofluorocarbons. 

Preferably the inventive preparations can also contain substances that absorb UV radiation in the 
UVB range, the total quantity of the filter substances being e.g. 0.1 wt.% to 30 preferably 
0.5 to 10 wt.%, especially 1 to 6 relative to the total weight of the preparation, in order to 
provide cosmetic preparations that protect the skin from the entire range of ultraviolet radiation. 
They can also act as sunscreens. 

The UVB filters can be oil-soluble or water-soluble. Oil-soluble substances that can be used 
include e.g.: 

- 3-benzylidene camphor derivatives, preferably 3-(4-methylbenzyliden) camphor, 3- 
benzylidene camphor; 

- 4-aminobenzoic acid derivatives, preferably 4-(dimethylamino)-benzoic acid (2-3thylhexyl) 
ester, 4-(dimethylamino)benzoic acid amyl ester; 

- Esters of cinnamic acid, preferably 4-methoxycinnamic acid (2-ethylhexyl) ester, 4- 
methoxycinnamic acid isopentyl ester; 

- Esters of salicylic acid, preferably 
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salicylic acid (2-ethylhexyl) ester, salicylic acid (4-isopropylbenzyl) ester, salicylic acid 
homomenthyl ester; 

- Derivatives of benzophenone, preferably 2-hydroxy-4-niethoxybenzophenone, 2-hydroxy-4- 
methoxy-4 ' -methy Ibenzophenone, 2,2 ' -dihydroxy-4-methoxybenzophenone ; 

- Esters of benzal malonic acid, preferably 4-methoxybenzal malonic acid die(2- 
ethylhexyl)ester; 

- 2,4,6-trianilino-(p-carbo-2'-ethyl-r-hexyloxy)-l,3,5-triazine. 

Water-soluble substances that can be used include e.g.: 

- Salts of 2-phenylbenzimidazole-5-sulfonic acid like its sodium, potassium, or triethanol 
ammonium salt, and sulfonic acid itself; 

- Sulfonic acid derivatives of benzophenones, preferably 2-hydroxy-4-methoxybenzophenone- 
5-sulfonic acid and its salts; 

- Sulfonic acid derivatives of 3-benzylidene camphor, such as e.g. 4-(2-oxo-3-bomylidene 
methyl) benzene sulfonic acid, 2-methyl-5-(2-oxo-3-bomylidene methyl) sulfonic acid and 
its salts. 

The subject-matter of the invention is also the combination of inventive active substances with 
one or a plurality of UVB filters and inventive cosmetic or dermatological preparations that also 
contain one or a plurality of UVB filters. 
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It can also be advantageous to combine the active substances with UVA fihers that were 
previously normally contained in cosmetic and/or dermatological. These substances are 
preferably derivatives of dibenzoylmethane, in particular l-(4'-tert. butylphenyl)-3-(4'- 
methoxyphenyl)propane-l,3-dion and l-phenyl-3-(4'-isopropylphenyl)propane-l,3-dion. These 
combinations and preparations that contain these combinations are also the subject-matter of the 
invention. The quantities used for the UVB combination can be employed. 

Advantageous preparations are obtained when the inventive active substances are combined with 
UVA and UVB filters. 

Combinations of the inventive active substances with one or a plurality of antioxidants and one 
or a plurality of UVA filters and/or one or a plurality of UVB filters are also particularly 
advantageous in accordance with the invention. 

The cosmetic or dermatological preparations can also contain inorganic pigments that are 
normally used in cosmetics for protecting the skin fi*om UV radiation. These are oxides of 
titanium, zinc, iron, zirconium, silicon, manganese, aluminum, cerium, and mixtures thereof, as 
well as variants in which the oxides are the active agents. Titanium dioxide-based pigments are 
particularly preferred. 

Subject-matter of the invention is also the method for producing the inventive topical 
preparation, which method is characterized in that the active substances are worked into cosmetic 
or dermatological formulations 
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in a manner known per se. 

Provided there is no information to the contrary, all data regarding quantities, proportions, and 
percentages relate to the weight and total quantity or to the total weight of the preparations. 

The following examples are intended to clarify the present invention without restricting the 
invention. 

Example 1 

Sun gel LF 4 (transparent) 

L-arginine hydrochloride 
Benzophenone-4 

Phenylbenzimidazole sulfonic acid 
Acrylamide/sodium acrylate copolymer 
Ethanol 
Glycerol 
NaOH(15%) 
Fragrance, preservatives 
Water, demin, (fully desalinated) 

Example 2 



Hydrodispersion 

Wt.% 

L-omithine hydrochloride 5.0 

Phenyl trimethicone 1 .0 

Carbomer (Carbopol 981) 1 .0 

Hydroxy propylmethyl cellulose 2.0 

Butylene glycol 3.0 

Tromethamine as necessary 



Wt.% 

10 

0.5 

1.3 

1.6 

5.0 

15.0 

as necessary 
as necessary 
to 100.0 
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EDTA solution (14%) 0.5 

Ethanol 5.0 

Fragrance, preservatives as needed 

Water, demin. (fully desalinated) to 100.0 

Example 3 

Suntan lotion OAV 

Wt.% 

L-arginine hydrochloride 5.0 

Octylmethoxycinnamate 5 .0 

Butyl methoxydibenzoylmethane 1 .0 

Cetearyl alcohol + PEG-40 castor oil + sodium cetearyl sulfate 2.5 

Glyceryl lanolate 1.0 

Lauryl methicone copolyol 0.5 

Mineral oil (DAB 9) 5.0 

Caprylic/capric triglycerides 5.0 

Acrylamide/sodium acrylate copolymer 0.3 

Cyclomethicon 2.0 

Ti02 1.0 

Glycerol 3.0 

EDTA solution (14%) 0.5 

Ethanol 5.0 

Fragrance, preservatives as necessary 

Water, demin. (fully desalinated) to 100.0 

Example 4 

Skin lotion W/0 

Wt.% 

L-arginine hydrochloride 2.5 

Folic acid 0.1 

Cyclomethicon 3.0 

PEG- 1 -glycerol sorbitan oleostearate 1.7 
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PEG-7 hvdrated ca^itor oil 




6.3 


Mineral nil /'HAR Q"^ 




13 0 


Caprylic/capric triglycerides 




13.0 


rrlvpprnl 

vJ 1 Y^^i \Ji 




4.0 


MgS04 




0.7 


r ragrance, preservatives 




do JlCL'Coodiy 


Water Hemin 


nxdinpic J 


to 100 0 






Wt.% 


T -nrcrininp VivdmplilnriHe 




5.0 






1.0 


ri-jKJ-J glyceryl blCmalC 




2 00 






3.00 


L^yciomeinicon 




1 00 


l^dpiy llv/Capi IL.* U IgljVx'Cl lUwi 




3.00 


l-'Ciyi dlCUIlUl 




3.00 


Oftvl mpthoxvrJnnamate 

\^\^\,y I iii^iiiw/\.Y villi lui licit V 




2.50 


EfUlallUl 




1.00 


Hyaluronic acid 




0.05 


1 ocopneryi dLeiaie 




0 50 


Glycerol 




4.00 


Fragrance, preservatives 




do UCWCooaiy 


waier, aemin. 


Examnle 6 


to 1 00 0 


w/u creme 




Wt.% 


L-citrulline 




2.5 


L-arginine hydrochloride 




2.5 


FAD 




0.1 


FMN 




0.05 
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PEG-22-dodecyl glycol copolymer 3.0 

Cetyl dimethicone copolyol 2.0 

Cyclomethicon 4.0 

Mineral oil (DAB 9) 4.0 

Caprylic/caprlc triglycerides 4.0 

Glycerol 4.00 

Fragrance, preservatives as necessary 

Water, demin. to 100.0 

Example 7 

W/0 creme 

Wt.% 

L-omithine hydrochloride 2.5 

L-citrulline 2.5 

L-arginine hydrochloride 2.5 

FAD 0.1 

FMN 0.05 

PEG-22-dodecyl glycol copolymer 3.0 

Cetyl dimethicone copolyol 2.0 

Cyclomethicon 4.0 

Mineral oil (DAB 9) 4.0 

Caprylic/capric triglycerides 4.0 

Glycerol 4.00 

Fragrance, preservatives as necessary 

Water, demin. to 100.0 

Example 8 

After-sun lotion 

Wt,% 

L-arginine hydrochloride 5.0 

Folic acid 0.1 

Urea 2.5 
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Cetearyl alcohol + PEG-40 castor oil + sodium cetearyl sulfate 2.50 

Glyceryl stearate SE 0.60 

Mineral oil (DAB 9) 4.00 

Caprylic/capric triglycerides 2.00 

Shea butter 2.00 

Avocado oil 2.00 

Tocopheryl acetate 3.00 

Acrylamide/sodium acrylate copolymer 0.30 

Glycerol 4.00 

Hyaluronic acid 0.05 

Bisabolol 0.05 

Fragrance, preservatives as necessary 

Water, demin. to 100.0 

Example 9 

Shower lotion 

Wt.% 

L-arginine hydrochloride 1 0 .0 

Folic acid 1.0 

Sodium laureth sulfate 12 

Cocamidopropyl betaine 5 

Cocamide DEA 1 

PEG-8 1 

Soybean oil 1 

Citric acid 0.1 

Sodium chloride 0.2 

Fragrance 0.1 

Water, demin. to 100.0 
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Example 10 

W/0 creme 





W I. /o 


L-orniinine nyurocnioriQc 


Z.J 


L-citrulline 




L^-arginine nyurocnioriae 


Z.J 


T7 A n 
r AL/ 


n 1 

U. i 


r iviiN 


U.UJ 


Folic acid 


n 1 


1 ocopneryi acciaie 


1 0 


urca 


^. J 


riiia-zz-aouecyi giycoi copoiynier 




Cetyl dimethicone copolyol 


2.0 


Cyclomethicon 


4.0 


Mineral oil (DAB 9) 


4.0 


Caprylic/capric triglycerides 


4.0 


Glycerol 


4.00 


Fragrance, preservatives 


as necessary 


Water, demin. 


to 100.0 
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Test report 



The excellent effect of the inventive active substances is demonstrated with the following test 
report. The lotion from Example 4, but with an active substance content of 5 w1:.% L-arginine 
and without folic acid and fragrance, is used in a "stinging test" with 5 wt.% lactic acid in water 
as "stinger", as in Frosch and Kligman (Peter J. Frosch and Albert M. Kligman: "A method for 
appraising the stinging capacity of topically applied substances", J. Soc. Cosmet. Chem. 28, 
197209, 1977). 

After one week's contra-lateral product application (two times daily), the stinging potential of 
5% lactic acid is determined after a single application to the cheeks. First on the one cheek and 
approx. 1 hour later, after the irritation has cleared, on the other cheek. The stinging potential of 
5% lactic acid is significantly reduced on the cheeks previously treated with 5% L-arginine. 
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Patent claims 

1 . Use of one or a plurality of compounds selected from the group of L-arginine, L-omithine, 
and L-citruUine or their salts, acid addition salts, esters, or amides, where necessary with the 
addition of folic acid or its salts and/or one or a plurality of compounds selected from the 
group of flavins, for prophylactic purposes and/or treatment of neurosensory phenomena. 

2. Use of cosmetic and dermatological topical preparations having a content of one or a 
plurality of compounds selected from the group of L-arginine, L-omithine, and L-citrulline or 
their salts, acid addition saUs, esters, or amides, where necessary with the addition of folic 
acid or its salts and/or one or a plurality of compounds selected from the group of flavins, for 
prophylactic purposes and/or treatment of neurosensory phenomena. 

3. Use in accordance with claim 2, characterized in that the preparations contain L-arginine. 

4. Use in accordance with claim 2, characterized in that the preparations contain folic acid 
and/or one or a plurality of compounds selected from the group of flavins. 
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5. Use in accordance with claim 2, characterized in that the preparations contain at least one 
antioxidant. 

6. Use in accordance with claim 2, characterized in that the preparations contain urea. 

7. Use in accordance with claim 2, characterized in that the preparations contain at least one 
UVA filter and/or at least one UVB filter and/or at least one inorganic pigment. 

8. Use in accordance with claim 2, characterized in that the preparations contain L-arginine, 
folic acid, one or a plurality of compounds selected from the group of flavins, at least one 
antioxidant, urea, and where necessary at least one UVA filter and/or at least one UVB filter 
and/or at least on inorganic pigment. 



